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LOCATING  THE  SOURCE  OF  LONG-DISTANCE  RADIO 
DISTRESS  SIGNALS  FROM  THE  SOUTHERN  OCEAN 


G.L.  Goodvln 


SUMMARY 


This  report  summarises  a  feasibility  study  and  strategy  appropriate  to 
search  and  rescue  operations  for  ships  or  aircraft  In  distress  in  the 
Southern  Ocean.  The  use  of  high  frequency  lonospherlcally-propagated 
radio  transmissions  Is  justified.  Measurement  of  the  direction  of 
arrival  of  a  radio  signal  facilitates  ray-tracing  through  the  ionosphere 
Characteristics  of  the  F2  region  mid-latitude  Ionisation  trough  are 
considered  In  detail.  Because  the  trough  Is  displaced  cowards  Australia, 
its  position  must  be  known  as  near  as  practicable  in  real  time  for 
sufficiently  accurate  ray-tracing  to  be  performed  quickly  enough  after 
receipt  of  a  distress  signal.  The  differential  Doppler  technique  using 
e.g.  150  and  400  MHz  transmissions  from  NAVSAT  (NNSS)  satellites  Is 
recommended  In  locating  the  trough  position. 
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1 .  FOREWORD 

The  resources  of  the  Australian  Defence  Forces  can  be  called  on  as  part  of 
a  search  and  rescue  strategy  In  an  expansive  region  adjacent  to  Australia, 

Including  the  Southern  Ocean  extending  south  and  southwest  from  Australia  to 
Antarctica.  The  present  report  summarises  a  feasibility  study  carried  out 
as  part  of  the  author's  study  programme  during  the  period  September,  1982  - 
March,  1983,  while  attached  to  the  Ionospheric  Studies  Group  at  the  Defence 
Research  Centre,  Salisbury  (DRCS)  on  leave  from  the  School  of  Physics  of  the 
South  Australian  Institute  of  Technology  (S.A.  Inst,  of  Tech.).  A  reasoned 
strategy  Is  presented  for  search  and  rescue  operations,  together  with  a 
condensation  of  relevant  background  Information. 

This  report  was  prepared  In  the  context  of  a  possible  Joint  proposal  between 
DRCS  and  the  S.A.  Inst  of  Tech,  to  apply  radio  direction-finding  techniques 
and  Ionospheric  ray  tracing  computations  In  determining  the  location  (ship/plane) 
of  long  range  radio  distress  signals.  This  proposal  did  not  proceed  because 
of  changing  priorities  within  the  Defence  Science  and  Technology  Organisation 
(DSTO) .  However,  Single  Station  Location  (SSL)  of  HF  distress  signals  propagating 
at  latitudes  containing  the  F2  region  midlatitude  trough  continues  to  pose 
a  real  problem,  the  nature  of  which  Is  examined  In  this  report. 

2.  ABSTRACT  OF  CONTENTS 

The  study  commenced  with  the  premise  Chat  ships/ 'planes  In  distress  will  have 
available  the  normal  means  of  radio  communication. 

It  Is  concluded  In  Section  3  that  over  the  next  several  years,  at  least,  a 
search  and  rescue  strategy  should  be  based  on  the  reception  In  Australia  of 
long  distance  lonospherically-propagated  transmissions  at  high  frequency  (H.F.) 
International  distress  frequencies  (l.e.  one  or  more  of  2.182,  4,125  and 
6.2135  MHz)  since  most  ships/ 'planes  have  comparatively  high  power  transmitters 
at  these  frequencies.  An  alternative,  such  as  the  use  of  low  power  emergency 
position  Indicating  radio  beacons  (EPIRB) ,  e.g.  at  121.5  MHz  to  be  received 
by  near-polar  orbiting  search  and  rescue  satellites  (SARSAT),  was  in  the 
proposal  stage  only,  when  this  report  was  compiled. 

In  Section  4,  it  is  assumed  that  the  direction  ot  arrival  (azimuth  and  elevation) 
of  lonospherically-propagated  signals  would  be  measured  by  means  of  at  least 
one  multi-element  receiving  aerial  array  In  southern  Australia.  An  extension 
of  the  existing  DRCS  RADFIN  computer  program  would  facilitate  ray-tracing  through 
the  Ionosphere,  In  computing  the  location  of  the  signal  source.  Section  4 
is  mainly  concerned  with  factors  that  limit  the  accuracy  with  which  a  signal  source 
can  be  located.  A  limitation  arises  through  a  lack  of  directly  observed 
Ionospheric  data  over  a  wide  expanse  of  Southern  Ocean.  A  further  limitation 
occurs  because  the  south  Invariant  magnetic  pole  (see  Appendix  1)  Is  displaced 
some  16  deg  of  latitude  from  the  geographic  pole  towards  Australia.  This  means 
that  phenomena  normally  expected  at  high  latitudes,  occur  nearer  to  Australia; 

In  particular,  the  F2  region  midlatitude  trough  of  Ionisation,  a  gross  Ionospheric 
feature,  can  occur  at  latitudes  as  far  equatorvard  as  the  Australian  continent. 
Although  the  trough  can  be  modelled  approximately,  as  considered  In  Section  4.3, 
its  width,  which  can  extend  through  10  deg  or  more  of  invariant  latitude  (defined 
In  Appendix  1),  is  unknown  without  direct  measurement.  In  addition,  the  trough's 
mean  position  is  a  function  of  magnetic  activity,  and  moves  equatorwards  under 
magnetically  disturbed  conditions;  unfortunately,  the  trough's  position  cannot 
readily  be  estimated  because  the  planetary  K  Index  (Kp) ,  which  Is  employed  as 
an  appropriate  measure  of  magnetic  activity.  Is  not  readily  available  until  long 
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after  required  for  computations.  It  is  concluded  that  the  accuracy  of 
ray-tracing,  and  hence  of  signal  source  position  location,  will  be 
severely  restricted  unless  the  midlatitude  trough  position  and  extent 
CcUi  be  measured  in  approximately  real  time.  The  use  of  the  Navy  Naviga¬ 
tion  Satellite  System  (NNSS)  and  a  differential  Doppler  technique  (Section 
7)  are  recommended  to  determine  the  trough  position  cir,d  latitudinal  extent. 

Section  5  summarises  features  of  some  satellite  systems,  with  particular 
reference  to  their  capacity,  or  otherwise,  to  serve  in  locating  ships/ 
'pleuies  in  distress.  In  particular,  the  MARISAT,  proposed  SARSAT,  NNSS 
and  Global  Positioning  System  (GPS)  systems  of  satellites  are  considered. 

In  Section  6,  the  characteristics  of  the  midlatitude  trough  are  discussed 
in  detail.  Its  general  location  and  features  are  considered,  including 
the  fact  that  although  the  trough  is  predominantly  a  winter  night-time 
phenomenon,  it  can  occur  unpredictably  at  other  times.  The  trough  forms 
just  equatorward  of  the  auroral  oval.  As  discussed  in  Appendix  3,  an 
important  factor  in  the  trough's  formation  is  the  transport  of  excited 
Nj  equatorwards  from  the  auroral  region  by  thermospheric  winds  and  by 
convection  associated  with  magnetospherically  produced  electric  fields. 

The  excited  Nj  undergoes  reactions  which  remove  ionised  atomic  oxygen 
cind  electrons  from  the  F2  region. 

The  Doppler  measurement  of  range  is  described  in  Appendix  2,  providing 
bac)cground  for  Section  7  in  which  the  differential  Doppler  technique  is 
outlined.  By  means  of  this  technique,  changes  in  phase  difference  between 
coherent  150  and  400  MHz  transmissions  from  a  NNSS  satellite  are  directly 
interpretable  as  changes  in  the  ionospheric  electron  content  along  the 
ray  path  from  satellite  to  receiver.  This  pennits  gross  changes  in 
electron  content  to  be  detected  as  the  ray  path  encounters  the  trough 
edges,  which  allows  the  trough's  position  and  extent  to  be  determined. 

Appendix  4  includes  the  results  of  preliminary  coii5>utations  of  model 
midlatitude  troughs,  using  a  modified  RADFIN  program. 


3.  AUSTRALIAN  EMERGENCY  COMMUNICATION  PROCEDURES 
3.1  Introduction 

The  Department  of  Defence  may  be  called  on  as  part  of  the 
land/sea  rescue  strategy  for  a  considerable  area  including  the 
Australian  continent,  a  large  section  of  the  Indian  Ocean  westwards 
of  Australia  and  a  region  south  and  southwest  of  Australia  to  the 
pole.  The  southern  region  includes  a  vast  expanse  of  the  Southern 
Ocean  which  is  the  particular  concern  of  the  present  report. 

Figure  1  (Ref.  26;  1980)  shows  graphically  that,  following  an  alarm, 
there  is  an  urgent  need  for  land/sea  rescue  procedures  to  be  initiated 
without  delay,  because  the  probability  of  human  survival  decreases 
rapidly  with  time.  This  figure  presumably  refers  to  "average"  clima¬ 
tic  conditions;  under  icy  conditions  at  high  latitudes,  survival 
times  could  be  considerably  shorter. 

The  use  of  high  and  medium  frequency  distress  signals  should  be  effec¬ 
tive  in  areas  such  as  the  Indian  Ocean  where  there  are  many  coastal 
stations  in  Australia,  India  and  Africa,  and  numerous  stiips  on  well- 
used  shipping  lanes.  These  can  all  assist  in  applying  direction¬ 
finding  techniques  to  locate  the  source  of  a  distress  signal. 
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Figure 


1.  Probeibility  of  survival  (%)  is  shown  as  a  function 
of  recovery  time  (from  Refs.  26  and  60) . 
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In  contrast,  much  of  the  Southern  Ocean  c«u:ries  little  commercial 
shipping,  cind  is  also  remote  from  coastal  stations  available  for 
the  reception  of  long-reuige  distress  signals.  For  this  reason, 
searching  operations  in  the  Southern  Ocean  would  rely  heavily  on 
direction-finding  by  me^uls  of  long-distance  transmissions  in  a 
region  where  published  ionospheric  characteristics  are  suspect 
because  of  the  edssence  of  directly  observed  ionospheric  data.  In 
particular,  the  properties  of  the  mid-latitude  "trough”  of  F-region 
ionisation  (e.g.  Ref.  22)  discussed  in  Section  4.2  can  cause  gross 
errors  in  direction  of  aurrival  observations  in  Australia  of  radio 
signals  propagated  via  the  ionosphere  from  a  source  south  of  Australia. 
The  characteristics  of  the  trough  are  largely  unmeasured  at  medium 
to  high  latitudes  in  the  Southern  Hemisphere,  and  must  be  inferred 
from  (limited)  Northern  Hemisphere  data  (e.g.  Ref.  39).  Direct 
measurements  of  the  trough  position  in  the  Southern  Hemisphere  would 
therefore  be  desirable,  as  a  function  of  time  of  day,  season,  sunspot 
number  and  the  planetary  index  (Kp)  of  magnetic  activity. 

Considering  the  size  of  this  region,  there  are  comparatively  few 
coastal  stations  available.  For  both  shipping  and  aircraft,  distress 
signals  would  usually  be  long-range  transmissions  via  the  ionosphere 
at  the  international  distress  frequency,  2.182  MHz,  or  the  supple¬ 
mentary  frequencies,  4.125  and  6.2155  MHz.  Maximum  transmitted  power 
is  400  W,  but  25  W  is  a  typical  value  for  small  craft.  (The  following 
discussion  largely  disregards  trcuismitter/receivers  operating  at 
higher  frequencies,  since  these  are  primarily  for  short-range  commun¬ 
ication  and  would  not  normally  involve  the  ionosphere.) 

Ref.  40  indicates  emergency  procedures  for  Australian  shipping, 
that  are  also  appropriate  for  aircraft,  (A  corresponding  book  is 
not  available  specifically  for  aircraft.) 

Ships  in  distress  rely  heavily  on  other  ship  and  shore  stations 
receiving  distress  signals  at  one  or  more  frequencies  reserved 
for  this  purpose.  Station  operators  are  obliged  to  keep  a  record 
of  p2u:ticulars  of  all  distress  messages  transmitted  or  received. 

A  distress  message  should,  if  practicable,  include  the  position 
(geographic  latitude  and  longitude  in  degrees  and  minutes)  of  the 
craft  in  distress,  or  its  position  in  relation  to  a  well-known 
geographical  feature.  In  the  preser .  discussion  it  is  assumed  that 
the  ship/ 'plane  location  is  either  unknown  or  requires  confirmation, 
so  that  the  measures  recommended  in  Section  4.1  need  to  be  put  into 
operation. 

3 . 2  Medium  and  High  Frequencies 

Ship  stations  are  expected  to  keep  the  maximum  watch  practicable 
on  the  international  radiotelephone  distress  and  calling  frequency, 
2.182  MHz,  particularly  during  the  three-minute  "silence"  periods 
reserved  exclusively  for  distress  signals,  conmencing  on  each  hour 
and  half-hour.  The  high  frequencies  (HF) ,  4.125  and  6.2155  MHz, 
are  supplementary  to  2.182  MHz  for  distress  and  safety  purposes, 
and  have  a  schedule  of  monitoring  by  the  Australian  0"erseas  Tele¬ 
communications  Commission  (OTC)  shore  stations.  (All  such  stations 
keep  watch  on  6.2155  MHz  within  their  station  hours,  which  are  at 
least  0500-1900  hours  local  time,  while  four  stations,  of  which 
only  Perth  is  relevant  to  the  Southern  Ocean,  keep  continuous  watch 
during  their  hours  of  service.) 

Where  practicable,  the  distress  ("Mayday")  signal  is  preceded  for 
one  minute,  by  an  alarm  signal  of  two  modulation  frequencies  trans¬ 
mitted  alternately  at  2200  and  1300  Hz. 
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The  following  propagation  guide  indicates  the  practical  limitations 
in  communication  from  ship-board  at  medium  and  high  frequencies: 

During  daytime,  when  ionospheric  D-region  2d5Sorption  is  strong, 

2  MHz  communication  is  limited  to  approximately  80  to  160  km, 
depending  on  the  ship's  transmitter  power  and  antenna  efficiency; 

2  MHz  is  the  first  choice  of  frequency  for  night-time  communication 
via  the  ionosphere. 

The  normal  day-time  frequency  is  4  MHz  for  distances  of  the  order 
of  300  km,  and  for  night-time  coimnunication  when  2  MHz  is  unsatis¬ 
factory  due  to  attenuation  with  distance  or  "static"  noise. 

For  day-time  coimnunication,  6  MHz  is  used  when  long  distances  make 
4  MHz  unsatisfactory.  The  use  of  6  MHz  is  not  recommended  at  night 
on  occasions  when  lower  frequency  communication  is  possible,  because 
of  the  interference  potentially  caused  and  received  with  long  dist¬ 
ance  6  MHz  transmissions. 

Single  sideband  transmission  at  2,  4  and  6  MHz  is  employed  mainly 
by  smaller  ships,  i.g.  fishing  vessels  and  non-commercial  ("pleasure") 
yachts  at  maxxmum  allowable  transmitter  powers  of  400  W.  In  prac¬ 
tice,  small  ship  users  normally  purchase  tr^u^smitters  with  100  W 
peak  envelope  power  which  typically  provide  25  W  of  aerial  power; 
transmitted  power  must  exceed  15  W  to  be  acceptctble  to  the  Common¬ 
wealth  Department  of  Communications.  Small  ships  commonly  have  a 
"wire"  (inverted  L)  aerial  or  a  "whip"  (approximately  vertical) 
aerial  if  the  ship  has  no  mast.  Several  HF  bands,  nominally  at 
2,  4,  6,  8,  12,  16,  22,  25  and  27  MHz,  are  in  general  use;  some 
b2uids  are  for  Morse  (i.e.  amplitude  modulation)  and  others  involve 
single  sideband  transmissions,  with  peak  envelope  powers  of  400  W. 
Their  use  is  for  long-range  communication  from  both  ships  and  air¬ 
craft.  Big  aircraft  use  continuously  variable  transceivers,  whereas 
light  aircraft  often  use  10  to  20  fixed  frequencies. 

Most  big  ships  would  normally  operate  at  high  powers  in  the  medium 
frequency  (MF)  maritime  bemd  (410-512  kHz)  for  long-distance  communi¬ 
cation  in  telegraphy  (Morse  code)  at  a  frequency  which  was  returned 
to  the  earth  by  the  ionosphere  both  in  the  day-time  and  at  night. 

They  would  generally  use  500  kHz  (or  else  2.182  MHz)  for  transmitting 
their  own  distress  signals  in  telegraphy.  Such  ships,  as  well  as 
shore  stations,  would  have  direction-finding  equipment  for  identi¬ 
fying  the  directions  of  signals  at  these  international  distress 
frequencies,  from  other  ships  in  distress. 

In  emergency  operations,  MF  and  HF  transmitters  on  craft  in  distress 
are  used  as  beacons  for  long-range  transmissions  over  several  hundred 
km  to  searching  aircraft  equipped  with  direction- finding  facilities. 

3.3  Very  High  Frequencies 

The  very  high  frequency  (VHF)  mauritime  mobile  radiotelephone  service 
operates  with  frequency  modulation  on  channels  in  the  ranges  156.375- 
156.80  MHz  and  161.575-161.95  MHz  providing  a  widely-used  service  for 
port  authorities  and  commercial  users,  e.g.  fishing  boats,  tugs  and 
oil  rigs,  as  well  as  for  a  growing  non-commercial  usage,  e.g.  for 
yachts.  In  peurticular,  156.80  MHz  is  allocated  as  an  international 
distress  and  safety  frequency  and  156.375  MHz  as  a  supple-sntary 
safety  frequency.  Ships  at  sea  are  expected  to  keep  a  listening 
watch  on  156.80  MHz.  The  emticipated  working  range  is  low,  being 
a  maximum  of  24  to  40  )un  for  ship  to  shore  communication,  and  13 
to  19  km  for  ship  to  ship,  corresponding  to  a  maximum  ship  trans¬ 
mitter  power  of  25  W,  and  subst^mtially  less  actual  power  transmitted. 
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Large  comniercial  'planes  over  the  Australian  mainland,  including 
flights  to  Tasmania,  are  in  constant  contact  with  air  traffic  control 
at  one  or  other  major  airport,  so  that  an  aircraft's  location  would 
be  known  reasonably  pecisely  in  case  of  emergency.  These  aircraft 
use  VHP  at  118  to  135.9  MHz  at  comparatively  low  powers  of  the  order 
of  several  watts.  Long  wire  transmitting-receiving  aerials  are 
commonly  used  in  aircraft,  e.g.  from  cockpit  to  tail,  or,  in  light 
aircraft,  a  horizontal  V-snaped  aerial  attached  to  the  tail,  cockpit 
and  wing. 

Emergency  positions  indicating  radio  beacons  (EPIRB)  are  being  progres¬ 
sively  fitted,  in  a  semipermanent  fashion,  to  Australian  ships  and 
light  aircraft.  An  EPIRB  is  a  battery-operated  low-power  radio 
transmitter  which  sends  a  repetitive  sawtooth-modulated  emergency 
signal  on  both  the  international  civil  and  military  aircraft  distress 
frequencies  (121.5  cmd  243  MHz  respectively).  Maximum  battery  life 
is  about  48  hours.  "Hand  held"  EPIRBs  are  portable  devices  carried 
e.g.  in  a  brief-case, which  are  more  common  in  light  aircraft  and 
could  he  used,  e.g.  in  rubber  dinghies.  The  aerial  is  often  a  "whip" 
about  31  cm  long  (a  quarter-wave  at  243  MHz) .  Alternatively,  flota¬ 
tion  EPIRBs,  e.g.  attached  to  a  life-raft  by  a  tow-line,  are  carried 
by  ships  and  airc^'aft  on  trans-oceanic  flights. 

It  has  been  claimed  (Ref.  26)  that,  world-wide,  there  are  approxi¬ 
mately  190,000  'planes  and  2000  ships  equipped  with  121.5/243  MHz 
EPIRBs.  This  does  not,  however,  reflect  the  situation  for  Australian 
'planes  on  domestic  flights  which  would,  in  general,  not  bo  equipped 
with  EPIRBs;  the  exception  would  be  for  flights  over  some  prescribed 
desolate  areas  of  the  Australian  interior.  EPIRP-  -id  to 

searching  aircraft,  equipp  d  with  dotection-findi.'.j  f-ilities,  in 
seeking  marine  and  aviation  casualt’  :s  in  and  around  Australia. 

The  detectable  range  of  an  EPIRB  transmission  at  sea  is  approximately 
100  km,  depending  on  the  height  of  the  searching  aircraft  and  its 
receiver  sensitivity. 

A  future  possibility  is  the  use  of  proposed  SARSAT  sateDites,  as 
considered  in  Section  3.2,  to  receive  low-power  EPIRB  transndssions 
and  to  use  Doppler  techniques  in  locating  the  transmission  source; 
this  is,  however,  in  the  proposal  stage  only. 

3.4  International  Distress  Frequencies  -  Summary 

As  detailed  in  the  above  discussion,  international  distress  frequ¬ 
encies  for  "May  day"  calls  are  2.182  MHz  (at  400  W  maximum  power 
for  SSB  transmissions),  27.880  MHz  (at  4  W  power),  156.8  MHz  (at 
25  W  power),  as  well  as  EPIRB  at  121.5  and  243  MHz.  Supplementary 
distress  frequencies  are  4.125  and  6.2155  MHz  (at  400  W) ,  27.860  MHz 
(at  4  W)  and  156.375  MHz  (at  25  W) .  In  addition,  the  500  kHz  band 
is  commonly  employed  for  distress  signals. 

For  practical  purposes,  direction- finding  for  vessels  in  distress 
in  the  Southern  Ocean  and  remote  from  the  Australian  coastline, 
would  be  restricted  to  HF  (2.182,  4.125  and  6.2155  MHz)  because  of 
the  comparatively  high  powers  transmitted  at  these  frequencies  and 
the  fact  that  they  can  be  propagated  over  long  distances  via  the 
ionosphere . 
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3.5  Radio  Installations  in  Australian  Shipping 

Table  1,  which  gives  information  on  the  frequencies  of  transmitter- 
receivers  installed  in  Australian  ships,  is  a  summary  of  data  provided 
from  the  Adelaide  office  of  the  Radio  Frequency  Management  Division, 
Australian  Department  of  Communications .  The  third  column  indicates 
that  in  those  States  (Victoria,  South  Australia,  Western  Australia 
and  Tasmania)  with  coastlines  on  the  Southern  Ocean,  5086  individual 
ships  have  HF  installations  at  2.182  MHz  capcdjle  of  long-range  com¬ 
munications;  in  contrast  the  (generally  smaller)  ships  which  belong 
to  clubs  are  equipped  predomincmtly  with  27  MHz  equipment  suitable 
for  limited  range  coast  and  ship  communication.  It  is  noted  that 
information  in  Table  1  is  restricted  to  commercial  and  "pleasure" 
craft  only. 

TABLE  1.  DISTRIBOTION  OF  COMMUNICATION  FREQUENCI.'S  IN  AUSTRALIAN  SHIPS 


INSTALLATION 

FREQUENCIES 

Ships  (generally  larger) 
equipp  l  with  radio¬ 
telephone  installations 
which  comply  with  the 
Wireless  Telegraphy  Act. 

2.182  MHz 

and 

possibly  also 
at  other 
frequencies 

5,086 

5,280 

27  MHz 

1,040 

2,118 

Ships  (generally 
smaller)  which  belong 
to  clubs,  e.g.  for 
pleasure,  fishing. 

27  MHz 

7,289 

HF,  VHF 
and  UHF 

92 

29 

Undefined 

2,206 

213 

4.  DIRECTION  OF  ARRIVAL  OF  RADIO  SIGNALS 

4.1  A  Facility  for  Locating  Sources  of  HF  Distress  Signals 

In  order  to  locate  the  sources  of  HF  radio  distress  signals  at 
middle  to  high  southern  invariant  geomagnetic  latitudes,  it  is 
proposed  that  at  least  one  multi-element  receiving  aerial  array 
be  established  in  southern  Australia,  e.g.  near  Adelaide,  to  act 
as  a  phase  interferometer.  The  aerial  outputs  of  the  array  would 
have  different  phases,  and  be  combined  to  produce  lobes  of  maximum 
reception  in  specified  directions.  This  would  allow  both  the  azimuth 
and  elevation  of  an  incoming  HF  signal  to  be  measured  to  a  specific 
accuracy.  One  or  more  arrays  of  this  )cind  are  expected  to  be  the 
basis  of  a  facility  whereby  sources  (ships/aircraft)  of  comparatively 
high  power  signals  can  be  located  rapidly  in  the  Southern  Ocean. 

Errors  can  arise  in  HF  direction-finding  because  of  the  occurrence 
of  ionospheric  irregularities  which  distort  the  otherwise  smoothly 
"reflecting'  ionosphere  and  cause  ionospherically-propagated  signals 
to  deviate  from  great-circle  paths  and/or  to  be  returned  from 
unexpected  heights.  Normally,  irregularities  arise  from  three  main 
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causes,  namely  travelling  ionospheric  disturbances,  ionospheric 
tilts  and  sporadic  E  (e.g.  Refs.  35  and  41).  However,  at  the  high 
invariant  geomagnetic  latitudes  of  the  Southern  Ocean,  the  mid-lati¬ 
tude  trough,  when  it  occurs,  will  produce  steep  gradients  of  electron 
density  which  could  be  the  major  source  of  errors  in  direction¬ 
finding  at  unpredictable  times. 

As  part  of  the  facility,  it  is  proposed  to  set  up  a  differential 
Doppler  satellite-receiving  system  to  obtain,  almost  in  real  time, 
the  location  of  the  F-region  mid-latitude  ionisation  "trough" 
(discussed  in  detail  in  Section  4.4  and  Appendix  4)  in  order  to 
facilitate  realistic  computations  of  the  ionospheric  ray  path  between 
source  and  receiver. 

In  summary,  the  proposed  facility  would  have  the  following  advantages: 

(a)  It  would  permit  rapid  determination  of  the  location  of  a  craft 
transmitting  HF  distress  signals  in  the  Southern  Ocean  where 
remoteness  from  major  shipping  lanes  and  the  lack  of  conven¬ 
iently  situated  shore-receiving  stations  present  severe  diffi¬ 
culties  in  locating  the  source  of  any  distress  signal. 

(b)  HF  transmissions,  in  particular  the  international  distress 
frequency,  2.  182  MHz,  and  the  supplementary  distress  frequen¬ 
cies,  4.125  cuid  6.2155  MHz,  are  available  from  virtually  all 
shipping  and  aircraft  with  long-distance  communication  facili¬ 
ties.  A  further  advantage  is  that  comparatively  high  power 
transmissions  are  employed,  up  to  a  maximum  of  400  W. 

(c)  In  the  well-established  Navy  Navigation  Satellite  System  (NNSS) , 
five  near-polar  orbiting  satellites  provide  very  high  frequency 
(VKF)  transmissions  at  150  and  400  MHz.  This  allows  a  set  of 
ionospheric  electron  content  measurements,  and  hence  an  update 
in  the  mid-latitude  trough  position,  to  be  made  at  comparatively 
short  mean  time  intervals  of  82  min  at  a  35  deg  latitude  ground¬ 
receiving  station. 

(d)  Use  of  NNSS  satellites  would  have  a  further  "spin-off"  advantage 
in  facilitating  ionospheric  research  over  the  Southern  Ocean. 

This  is  a  region  which  is  interesting  because  of  its  proximity 
to  Australia,  is  complicated  by  the  mid-latitude  trough,  and 
where  the  ionosphere  has  hardly  been  studied  previously  due 

to  a  lack  of  islands  on  which  to  establish  land-based  ionospheric 
soui.ding  stations. 

The  highest  latitude  at  which  a  trough  might  be  observed  in  Adelaide 
using  the  differential  Doppler  technique  is  indicated  from  the  follow¬ 
ing  arguments.  It  is  assumed  that  a  representative  trough  has  a 
peak  elevation  of  about  450  km  (justified  in  this  sub-section  below) 
and  has  an  apprecieUble  effect  on  total  electron  content  measured 
from  satellite  signals  passing  through  the  ionospheric  trough  at 
elevations  up  to  450  km.  (A  trough  with  a  rather  greater  peak  eleva¬ 
tion  of  600  km  is  illustrated  in  Figure  2  from  Ref.  11) .  It  is 
further  assumed  that  at  elevation  angles  of  10  deg  (ard  above)  signals 
are  received  in  Adelaide  with  strengths  large  enough  to  be  useful. 

From  geometry,  a  trough  which  was  just  observable  at  a  10  deg  elevation 
angle  would  therefore  occur  at  a  geographic  latitude  approximately 
13  deg  south  of  Adelaide  (lat.  35  S) ,  i.e.  at  a  geographic  latitude 
approximately  48  S. 
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As  discussed  in  Section  4.1,  the  mean  position  of  troughs  is  aibout 
60  deg  invariant  latitude  which,  for  Adelaide's  longitude,  corres¬ 
ponds  to  a  geographic  latitude  of  approximately  48  S;  (see  Figure  3). 
It  follows  that  troughs  which  occurred  approximately  in  the  mean 
trough  position  would  be  observ2d}le  in  Adelaide.  Furthermore, 
because  of  their  appreciable  latitudinal  extent,  troughs  which  were 
centred  at  higher  latitudes  but  which  extended  through  the  48  S 
invariant  latitude  should  also  be  observable  in  Adelaide.  Troughs 
should  cert^tinly  be  observable  when  they  were  located  sufficiently 
equatorwards  of  the  mean  position  so  as  to  produce  maximal  effects 
on  ionospherically-propagated  transmissions  from  the  Southern  Ocean, 
which  is  the  point  of  main  concern  in  the  present  context. 

In  attempting  to  locate  a  craft  in  distress  as  rapidly  as  possible, 
the  following  procedure  is  proposed: 

(a)  Information  that  an  HF  distress  signal  has  been  received  is 
forwarded,  e.g.  by  OTC,  without  delay  to  the  person  responsible 
for  operating  the  position-finding  facility. 

(b)  The  direction  of  arrival  of  the  radio  signal  is  measured  in 
both  azimuth  and  elevation,  at  a  receiving  station  near  Adelaide 
(and  possibly  elsewhere)  using  a  highly  directional  multi¬ 
element  crossed  aerial  array. 

(c)  A  modification  of  the  DRCS  RADFIN  ray-tracing  program  is 
employed  to  compute  an  estimate  of  the  signal  source  position, 
which  is  to  be  relayed  to  the  search  and  rescue  authority. 
Currently,  RADFIN  computations  are  based  on  CCIR  data  in  Ref.  14 
constituting  a  representative  ionosphere.  These  data  require 
modification  and  improvement  to  form  a  composite  model  iono¬ 
sphere,  by  superimposing  on  the  CCIR  data  the  best  available 
definitive  model  of  the  mid-latitude  trough,  as  devised  by 
Halcrow  and  Nisbet  (Ref.  39).  in  their  model,  the  trough  posi¬ 
tion  is  predicted  in  invariant  geomagnetic  co-ordinates  (Appen¬ 
dix  1  and  Ref.  7)  by  utilising  the  parameters,  universal  time, 
day  of  year  (i.e.  season),  sunspot  number  and  planetary  K  index 
(Kp)  which  is  a  measure  of  world-wide  geomagnetic  activity 
(Ref.  47)  . 

In  practice,  Kp  cannot  reliably  be  obtained  sooner  than  some 
four  to  six  weeks  after  observation.  In  addition,  Halcrow  and 
Nisbet' s  model  is  based  on  ionospheric  observations  in  the 
Northern  Hemisphere  where  the  north  invariant  geomagnetic  pole 
is  at  approximately  80  deg  N  geographic  latitude,  in  contrast 
to  the  south  invariant  pole  at  approximately  74  deg  S  geographic 
latitude  (Ref.  28) .  The  composite  model  ionosphere  would  be 
based  on  only  an  estimate  of  Kp.  It  can  therefore  be  regarded 
only  as  a  first  approximation,  but  it  is  useful  to  insert  into 
the  program  to  obtain  an  initial  estimate  of  the  location  of 
the  signal  source.  This  information  is  relayed  to  the  search 
and  rescue  authority. 

(d)  Although  the  effect  of  the  mid-latitude  trough  is  greatest  on 
winter  nights,  its  effect  at  other  times  is  also  appreciable 
and  unpredictable  (Ref.  15).  As  a  next  step,  direct  observa¬ 
tions  of  the  location  of  the  trough  are  therefore  made  by  means 
of  NNSS  satellites  in  employing  the  differential  Doppler  effect. 
This  permits  a  direct  measurement  of  the  trough  position  and 
latitudinal  extent  to  be  inserted  into  the  progreun,  to  replace 
the  estimated  trough  position  previously  computed  within  the 
program,  from  parameters  w)\ich  included  an  unreliable  Kp  value. 
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Figure  3.  Invariant  Geomagnetic  Co-ordinates,  Epoch  1969.75 
(from  Evans  et  al  (1969),  Ref.  28)  . 
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Furthermore,  the  variation  in  total  electron  content  observed 
during  the  passage  of  a  satellite  could  be  translated  into 
information  on  the  vertical  profile  of  ionospheric  electron 
density  in  the  trough.  This  should  permit,  for  the  Southern 
Hemisphere,  the  eventual  improvement  of  the  trough  model 
currently  available  for  the  program  (Ref.  39)  which  is  based 
on  Northern  Hemisphere  observations. 

(e)  The  program,  with  the  measured  trough  position  included,  is 

then  employed  to  compute,  more  accurately,  the  location  of  the 
signal  source;  the  information  is  relayed  to  the  search  and 
rescue  authority. 

4.2  Maximal  Direction  of  Arrival  Errors 

The  following  comments  are  mainly  restricted  to  propagation  problems 
associated  with  the  mid-latitude  trough.  It  is  noted  that  problems, 
particularly  related  to  direction  of  arrival  measurements,  disappear 
when  propagation  is  via  a  sporadic  E  (Eg)  layer,  from  which  reflec¬ 
tions  would  be  predominantly  specular. 

Buchau  et  al  (Ref.  12)  used  Alouette  topside  ionogram  data  as  a 
basis  for  three-dimensional  ray-tracing  of  radio-waves  assumed  to 
originate  from  the  Polar  Fox  II  over-the-horizon  radar  system,  at 
the  site  indicated  as  a  source  of  “rays”  in  Figure  4.  The  isoionic 
contours  in  this  figure  indicate  the  location  of  the  (cross-hatched) 
mid-latitude  trough.  The  figure  gives  some  indication  of  strong 
north-south  electron  density  gradients;  ground-based  and  airborne 
measurements  show  that  steepest  gradients  are,  in  fact,  greater 
than  the  figure  would  suggest,  being  approximately  2  to  3  MHz  per 
100  loll. 

Figure  5  from  Ref.  12  shows  a  set  of  curves,  each  of  which  is  the 
locus  of  landing  points  of  5  MHz  radio-waves  transmitted  in  a  parti¬ 
cular  azimuthal  direction  at  various  elevations.  Zero  azimuth 
corresponds  to  transmissions  in  a  direction  45  deg  E  of  N  which  is 
approximately  along  the  trough.  Figure  5  indicates,  for  instance, 
that  rays  directed  at  approximately  40  deg  E  of  N  (i.e.  at  -5  deg 
to  the  length  of  the  trough)  which  land  at  ship  distances  of  the 
order  of  3000  km  are  some  1400  km  or  28  deg  off  a  great  circle  path 
as  a  result  of  the  presence  of  the  mid-latitude  trough.  This  is 
ted^en  to  be  a  measure  of  the  maximal  error  that  might  arise  in 
attempting  to  locate  a  source  of  transmissions  by  measuring  the 
direction  of  arrival  of  the  signal,  unless  allowance  is  made  for 
the  presence  of  the  mid-latitude  trough. 

4.3  Limitations  in  Ionospheric  Data  for  the  Southern  CX:ean 

Ionospheric  ray-tracing  techniques  at  HF,  with  particular  applica¬ 
tion  to  travelling  ionospheric  disturbances,  are  discussed  in 
Ref.  35. 

By  means  of  the  RADFIN  computer  progreun  developed  at  ORCS,  the 
location  of  a  source  of  transmissions  received  via  the  ionosphere 
may  be  computed  from  measurements  of  the  direction  of  arrival  of 
signals.  Ray-tracing  techniques  are  applied  to  ionospheric  data, 
within  the  limits  of  the  most  con^rehensive  data  set  available 
(Ref.  15).  The  extent,  and  limitations,  of  these  CCIR  data,  wtU.ch 
are  employed  in  the  calculation  of  world-wide  F2  region  ionospheric 
characteristics,  are  indicated  by  the  following  description  from 
Ref.  15: 


00 


Figure  4.  Contours  of  Median  ^^2  (in  MHz)  Derived  from 
Alouette  1  Topside  Soundings.  Data  base: 
Winter  1963,  0000  to  1200  UT,  Kp  0  to  1+.  The 
co-ordinate  system  is  corrected  geomagnetic 
latitude  and  corrected  geomagnetic  local  time 
(from  Buchau  et  al  (1973) ,  Ref.  12)  . 
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Ray tracing  Results  for  5  MHz.  Shown  are  the  loci 
of  the  landing  points  as  a  function  of  the  initial 
azimuths.  Initial  azimuths  with  respect  to  a 
central  boresight  of  45  deg  east  of  north  arc 
indicated  by  the  nunijers  to  the  right  of  the  curves 
(from  Buchau  et  al  (1973),  Ref.  12). 
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"Data  for  five  years,  1954-  1958,  were  analysed  month  by  month. 

These  years  were  typical  of  a  period  from  solar  cycle  minimum  to 
maximum.  In  addition  to  data  from  stations  actually  operating 
in  those  years,  it  was  possible  to  use,  as  supplemental  data, 
observations  from  stations  that  were  in  operation  in  other  years, 
to  fill  some  large  gaps  in  world  coverage.  This  is  particularly 
helpful  for  1954,  a  period  of  solar  cycle  minimum.  The  use  of 
data  from  1964,  also  a  solar  cycle  minimum  year,  to  supplement 
the  1954  data  was  necessary  to  provide  sufficient  data  for  mapping 
most  of  the  Southern  Hemisphere.  The  severe  limitations  in  Southern 
Hemisphere  data  are  particularly  noted.” 

Although  the  CCIR  data  indicate  a  smoothed  "average"  position  of 
the  trough,  as  seen  in  Figures  6  and  7,  they  do  not  include  steep 
ionisation  gradients  that  are  important  in  computing  directions  of 
arrival  of  signals  refracted  by  the  ionosphere. 

Limitations  in  the  accuracy  of  CCIR  numerical  maps  are  indicated 
by  the  following  statement  from  Ref.  14.  Comparisons  made  (Refs. 

13  and  44)  "between  satellite  observations  of  plasma  frequencies 
and  the  numerical  maps  of  foF2  reveal  that  average  discrepancies 
are  generally  small,  but  that  in  some  ocean  regions,  notadbly  the 
South  Pacific  and  in  the  Southern  Hemisphere,  and  at  some  times, 
the  predictions  are  in  error  by  as  much  as  a  factor  of  2". 

In  Ref.  16  it  is  comnented  that  "the  ionosonde  network  in  existence 
when  Report  340  (CCIR  (1967))  was  prepared  has  provided  a  reasonable 
sample  of  one-third  of  the  Earth's  surface.  Since  that  time,  approxi¬ 
mately  15  new  ionosondes  have  begun  operation..." 

However,  the  major  source  of  error  in  ionospheric  data  for  the 
Southern  Ocean  still  arises  from  the  lack  of  ionospheric  sounding 
stations  in  this  region  south  and  south-west  of  Australia.  This 
necessitates  the  extrapolation  between  data  obtained  at  stations 
thousands  of  km  apart,  as  seen  in  Figure  8.  The  limited  accxiracy 
of  data  for  the  Southern  Ocean  is  highlighted  by  the  absence  of 
stations  in  cm  area  bounded  by  the  geographic  latitudes  of  Hobart 
(42.9  S)  and  Antarctica  (e.g.  Wilkes  66.9  S;  Mirny  66.5  S;  Terre 
Adelie  66.7  s)  cmd  the  longitudes  of  Macquarie  Island  (159.0  E) 
and  Kerguelen  (70.2  E)  (Ref.  21).  These  stations  are  listed  in 
T^lble  2  together  with  other  adjacent  stations. 

4.4  Usefulness  and  Limitations  of  the  Mid-latitude  Trough  Model 

In  order  to  make  meaningful  ray-tracing  computations,  a  realistic 
model  of  the  trough  must  be  superimposed  on  the  CCIR  data.  In  the 
present  context,  this  has  been  possible  by  adding  to  the  RADFIN 
program  subroutines  obtained  from  the  Pennsylvania  State  University 
(Mr.  R.L.  Divany)  which  are  based  on  the  mid-latitude  trough  model 
developed  by  Halcrow  and  Nisbet  (Ref.  39). 

In  the  present  discussion,  it  is  assumed  that  this  model,  derived 
for  the  Northern  Hemisphere,  can  be  applied  directly  to  the  Southern 
Hemisphere.  Ref.  78  provided  some  measure  of  support  for  this 
assumption  from  Ariel  III  satellite  observations  at  conjugate  regions 
in  the  two  hemispheres  by  showing  that  the  low  latitude  boundaries 
of  troughs  in  the  two  hemispheres  coincided  within  ±1.2  deg  of 
invariant  latitude  under  a  variety  of  conditions. 

Other  minor  improvements  to  the  model  were  recommended  by  Loc)cwood 
(Ref.  54);  these  are  briefly  considered  later  in  this  section. 


Figure  6.  Arrows  indicate  the  mid-latitude  trou^  position 
(from  CCIR  REPORT  340  (1966),  Ref.  15). 
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iqure  7.  Arrows  indicate  the  mid-latitude  trough  position 
(from  CCIR  REPORT  340  (1966),  Ref.  15). 


southern  Ocean  adjacent  to  Australia 
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TABLE  2.  LOCATIONS  OF  HIGH-LATITUDE  SOUTHERN  STATIONS 
(Data  from  Conkright  and  Brophy  (1982J ,  Ret.  21) 


LAT. 
deo  s 

LONG, 
deq  E 

RELEVANT  lONOGRAM 
RECORDING  PERIOD  FOR 
CCIR  (1966)  DATA 

Campbell  Island 

52.50 

169.20 

1954-58,  1964 

Stations 

(Casey) 

(66.28) 

(110.53) 

(2/1969-1975) 

between 

Geographic 

Hobart 

42.90 

147.30 

1954-59,  1964 

Latitudes 

Kerguelen 

49.35 

70.24 

1954-58,  1964 

40  and  70  S 
and 

Macquarie  Island 

54.50 

159.00 

1954-11/1958 

longitudes 

Mirny 

66.50 

93.00 

4/1956-58,  1964 

70  and  170  E 

Terre  Adelie 

66.66 

140.02 

1954-58,  1964 

WiUces 

_  _ 

66.90 

110.50 

7/1957-58,  1964 

Canberra 

35.30 

149.00 

1954-58,  1964 

Cape  Hallett 

72.30 

170.20 

1957-58,  2/1964-12/64 

Christchurch 

43.60 

172.80 

1954-58,  1964 

Little  America 

78.20 

197.80 

7/1957-58 

Other 

Mawson 

67.60 

62.90 

2/1958,  1964 

(Adjacent) 

Mundaring 

31,99 

116.31 

1964 

Stations 

Salisbury 

34.70 

138.60 

1964 

Scott  Base 

77.90 

166 . 80 

3/1957-58,  1964 

Vostok 

78.4 

106.9 

3/1958-12/58,  1964 

Watheroo 

30.30 

115.90 

1954-58 

Woomera 

31.00 

138.6 

1964 
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The  position  of  the  mid-latitude  trough  shows  substantial  temporal 
variations  in  latitude.  Figures  9a  and  9b  are  examples  of  true 
height  versus  latitude  distributions  (derived  by  G.D.  Slade  and 
I.P.  Buttery  from  CCIR  data  using  the  RADFIN  program).  The  trough 
appears  in  these  figures  as  only  a  small  change  in  height  over  a 
broad  range  of  latitudes  and  centred  roughly  on  the  mean  latitude 
expected  for  the  trough. 

To  represent  the  iono  jhere  in  the  absence  of  the  trough,  Halcrow 
(Ref.  38)  and  Halcrow  and  Nisbet  (Ref.  39)  employed  F2  layer  pea)c 
electron  densities  (denoted  Nj,F2(CCIR))  derivable  from  tabulated 
CCIR  data  based  on  observations  (Ref.  15) . 

In  effect,  this  necessitated  interpolating  between  data  representing 
sunspot  numbers  0  and  100.  It  is  noteworthy  that  this  procedure 
also  introduced  a  potential  source  of  error  by  neglecting  the  fact 
that  the  CCIR  data  already  included  the  influence  of  the  mid-latitude 
trough  averaged  over  all  latitudes  where  it  occurred.  The  trough 
occurred,  in  fact,  at  lower  latitudes  for  higher  sunspot  numbers, 
corresponding  to  enhanced  geomagnetic  activity.  A  comparison  of 
the  particular  CCIR  data  in  Figures  6  and  7  shows  that  the  geogra¬ 
phic  latitude  of  the  average  trough  position,  south  of  Australia, 
for  sunspot  number  zero  was  some  6  deg  polewards  of  its  average 
position  for  sunspot  number  100.  Figures  9a  and  9b,  which  show  true 
height  as  a  function  of  geographic  latitude,  also  indicate  that  the 
trough  moved  equctorward  by  several  degrees  of  latitude  as  the  sun¬ 
spot  number  increased  from  0  to  100. 

Halcrow  and  Nisbet  (Ref.  38)  indicated  that  in  their  model,  in  the 
particular  case  of  the  spring  equinox,  the  northern  and  southern 
bottom  edges  of  the  trough  represent  peak  f2  region  electron  densi¬ 
ties  (N[nF2)  wliich  ate  0.25  of  the  densities  in  the  absence  of  the 
trough,  as  seen  in  Figure  10.  They  employed  29770  Northern  Hemis¬ 
phere  observations  of  peak  F2  layer  electron  densities,  denoted 
N^F2  (observed),  obtained  by  means  of  Alouette  I  and  II  (Refs.  4 
and  5)  satellite  topside  sounders  from  1962  to  1970.  They  used  the 
ratio: 

4)  =  Nnf  2  (observed) /iN4^2  (CCIR) 

at  the  intersection  of  invariant  latitude  lines  and  the  300  km  alti¬ 
tude  level,  in  defining  the  locations  of  the  walls  of  observed 
troughs,  as  a  basis  for  developing  their  trough  rodel. 

In  Ref.  39,  Halcrow  and  Nisbet  stated  that  valu -s  of  ()>  at  the 
northern  and  southern  bottom  edges  were  similar  3  each  other,  and 
exhibited  no  significant  variation  with  season,  Kp,  sunspot  number 
or  local  time.  In  particular,  they  indicated  that  a  constant  value 
of  4>  =  0.25  ±  20%  provided  as  satisfactory  a  fit  to  the  available 
data  as  did  more  coit^licated  expressions. 

In  order  to  obtain  some  estimate  of  the  error  involved  in  applying 
Halcrow  and  Nisbet's  model  (Ref.  39),  the  particular  case  of  winter 
(July)  and  a  local  time  corresponding  to  20.00  O.T.  at  geographic 
longitude  225  M  are  considered.  It  is  assumed  that  the  model,  which 
was  derived  for  the  Northern  Hemisphere,  also  applies  to  the  Southern 
Hemisphere.  These  conditions  pertain  to  Figures  6  and  7  which  incxude, 
in  the  top  maps,  contours  of  the  F2  region  extraordinary  ray  critical 
frequency  (fjjF2)  at  vertical  incidence.  These  show,  in  particular, 
that  at  geographic  latitude  50  S,  fjjF2  =  2.6  MHz  in  both  figures, 
whereas  at  60  S,  =  2.4  MHz  for  sunspot  number  zero  and 

FjjF2  =  3.0  MHz  for  sunspot  number  100.  These  latter  values  of  fj(F2 
at  60  S  would  correspond  to  an  increase  in  electron  density  of 


U.T  Tim*  1M0  Sun  Spot  Ho,  1M 


Curves  of  constant  ionisation  density  computed  with  the  DRCS  "RADFIN" 
program,  by  G.D.  Slade  and  I.P.  Buttery  for  the  longitude  of  Adelaide,  July. 
Plasma  frequencies  (MHz)  of  curves  are  indicated. 


program,  by  G.D.  Slade  and  I.P.  Buttery  for  the  longitude  of  Adelaide, 
Plasma  frequencies  (MHz)  of  curves  are  indicated. 
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approximately  50%  as  sunspot  number  changed  from  0  to  100.  It  is 
concluded  that  veuriations  in  F2  region  peeik  electron  density  assoc¬ 
iated  with  sunspot  number  v2u:iations  can  potentially  cause  consider- 
cdDle  errors  in  trough  positions  computed  on  the  basis  of  a  constant 
value  of  (j)  =  0.25.  In  particular,  the  possible  error  in  (j)  of  ±20% 
claimed  in  Ref.  39,  might  well  be  an  optimistic  underestimate. 

Halcrow  and  Nisbet's  model  predicts  values  of  ij)  which  are  used  as 
multipliers  in  converting  values  of  peeOc  F2  region  electron  density 
observed  when  the  trough  is  present.  It  is  noted  that  this  model 
provides  values  corresponding  only  to  pecUt  electron  densities. 

In  ray-tracing,  e.g.  in  using  the  RADFIN  program,  it  is  necessary 
to  generate  an  electron  density  profile  representative  of  the  F2 
region  trough  extending  below  its  peak,  i.e.  values  of  iji  that  are 
appropriate  to  various  heights  below  tiie  peak  are  produced,  as 
discussed  below  in  this  section. 

Exeunples  of  trough  profiles  are  indicated  in  Figure  11  from  Ref.  12; 
the  corresponding  directions  of  radar  transmissions  can  be  inferred 
from  Figure  i2,  as  can  the  fact  that  the  second  and  third  profiles 
in  Figure  11  represent  rays  which  passed  twice  through  the  trough, 
la  Figure  11,  t.he  trough  which  has  a  critical  freguency  (foF2)  of 
approximately  2  MHz,  is  about  100  and  300  km  above  the  corresponding 
electron  density  contours  on  its  equatorward  and  poleward  sides, 
respectively.  The  trough  may  therefore  be  represented  as  having 
a  height  of  eibout  200  km.  Its  altitude  is  approximately  450  km  at 
its  pe20c. 

Figures  2,  13,  14  and  15  from  Ref.  11  show  that  the  height  of  the 
peak  electron  density  occurred,  on  four  occasions,  approximately 
450,  360,  450  and  350  km  respectively  above  the  level  of  a  similar 
electron  density  at  latitudes  adjacent  to  the  trough.  A  representa¬ 
tive  trough  height  of  about  400  km  therefore  seems  appropriate.  The 
trough  critical  frequencies  <f<jF2)  range  from  1.5  to  2.2  MHz,  and 
the  trough  altitudes  are  approximately  630  km  at  their  peaks. 

Data  from  the  Advanced  Ionospheric  Sounder  (AIS)  are  presented  in 
Ref.  63,  including  two  representative  troughs  shown  in  Figure  16. 
These  have  heights  of  approximately  360  euid  100  km  and  respective 
fQF2  values  of  2  and  1.5  MHz.  The  troughs  have  altitudes  of  approxi¬ 
mately  600  and  430  km  respectively  at  their  peaks.  Further  AIS  data 
are  presented  by  Dudeney  et  al  in  Ref.  23,  which  includes  Figure  17 
representing  a  trough  of  height  approximately  200  km,  f(jF2  equal  to 
2  MHz,  and  an  altitude  of  approximately  470  km  at  its  peak. 

Lockwood  (Ref.  53)  measured  the  Doppler  shifts  of  7.335  MHz  signals 
transmitted  over  great  circle  paths  of  5212  km  such  that  double  hop 
F2  echoes  could  be  returned  to  ground  from  the  trough  in  two  regions 
approximately  one-quarter  and  three-quarters  of  the  distance  from 
transmitter  to  receiver.  Using  a  siitple  ray-tracing  model,  and 
assuming  "reflection  points  at  the  low  altitudes  of  200  km",  he  was 
able  to  interpret  his  observations  in  terms  of  two-hop  F2  echoes. 
Hence,  he  deduced  that  the  height  of  the  F2  region  peak  electron 
density  within  the  trough  was  only  of  the  order  of  30-80  km  above 
the  height  having  this  electron  density  outside  the  trough.  This 
conclusion  was  in  accord  with  Lockwood’s  further  report  (Ref.  54) 
on  focusing  gains,  or  losses,  of  up  to  30  dB  for  double-hop  F2 
region  reflections  via  the  trough  over  these  long  paths. 


Continuous  Aurora  (shown  in  Fig.  4  of  Ref.  12)  Transferred 
into  the  Corrected  Geomagnetic  Co-ordinate  System.  Note 
extremely  good  representation  of  the  aurora  midnight 
sector  in  corrected  geomagnetic  co-ordinates.  Black 
square  indicates  location  of  POLAR  FOX  II.  The  sector 
of  the  ionosphere  studied  lies  between  -38  deg  and 
+38  deg  azimuth  with  respect  to  the  0  deg  boresight  in 
25  radial  cross-sections  (after  J.  Buchau  et  al  (1973), 
Ref.  12). 


ERL-0449-SD 


30 


STAUCTURC  OF  ThC  MMN  TIUUGH 


Figure  16.  l3o-ionic  Contours  of  Different  Types  of  Mid-latitude 
Trough.  Numbers  associated  with  contours  indicate 
plasma  frequency  in  MHz  (after  Rodger  and  Pinnock 
(1980),  Ref.  63). 
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It  is  noted  that  at  a  frequency  of  f  =  7.335  MHz  and  on  Lockwood's 
(Ref.  5  3)  assvm4>tion  of  a  200  km  reflection  altitude  for  two-hop 
echoes  and  a  5212  km  path,  the  equivalent  vertical  frequency  (fg^) 
(e.g.  Ref.  9)  which  equals  the  plasma  frequency  at  the  reflection 
height,  would  be,  neglecting  the  earth’s  curvature: 

fgy  =  fcosi 

where  i  =  cot”' (200/1303)  is  the  angle  of  incidence  of  the  signal 
on  the  ionosphere.  Hence,  fg^  -  1.1  MHz  wliich  being  below  the 
values  of  {^^2  expected  for  the  trough,  implies  that  the  signal 
would  be  returned  to  earth. 

From  the  foregoing  examples  of  troughs,  450  km  would  seem  to  be  a 
reasonable  estimate  for  the  altitude  of  the  peak.  On  the  assumption 
of  450  )cm  as  a  possible  reflection  altitude  for  two-hop  echoes,  it 
would  follow  that  i  =  cot"' (450/1303) ,  and  hence  f^^  =  2.4  MHz. 
However,  the  values  of  fQF2  for  the  troughs  ranged  from  1.5  to 
2.2  MHz,  so  that  7.335  MHz  signals  for  which  fgy  =  2.4  MHz  would, 
in  fact,  have  penetrated  these  troughs.  This  simple  argument  there¬ 
fore  casts  doubt  on  Lockwood's  (Refs.  53  and  54)  estimate  of  50  km 
for  trough  heights.  It  seems  at  least  possible  that  his  signals 
may  have  been  reflected  from  the  walls  of  troughs  rather  than  from 
their  peaks. 

On  the  basis  of  his  observations,  Lockwood  (Ref.  54)  made  the 
following  minor  modifications  to  Halcrow  and  Nisbet's  (Ref,  39) 
model ! 

(1)  To  avoid  the  effects  of  discontinuities  in  the  gradient  of 
at  the  top  and  base  of  each  wall,  the  linear  walls  were 
replaced  by  a  cosinusoidal  increase  up  to  the  constant  maximum 
value  of  $. 

(2)  Corresponding  to  the  trough  in  the  F2  region  peak  electron 
density,  a  similar  increase  was  included  in  the  height  of  the 
F2  region  peak  by  raising  it  by  an  amount  Ah  inside  the  trough. 
The  model  variation  of  Ah  was  made  the  same  as  for  i}>!  in  parti¬ 
cular,  across  the  trough  bottom.  Ah  assumed  a  maximum  value 

of  50  Ion  coincident  with  the  maximum  depletion  in  electron 
density  for  which  Halcrow  and  Nisbet  (Ref.  39)  assumed  that 
<()  =  0.25  ±  0.5. 

(3)  The  F2  layer  semi  thicknesses  were  assumed  to  be  constant  across 
the  trough . 

Rush  et  al  (Ref.  66)  have  published  a  simplified  model  of  the  high 
latitude  ionosphere  for  telecommunications  applications.  It  employs 
CCIR  data  to  predict  median  monthly  values  of  the  critical  frequency, 
height  of  the  maximum  electron  density  and  semi thicknesses  of  the 
E,  FI  euid  F2  regions.  It  also  predicts  high  latitude  features 
including  the  F2  region  ionisation  trough,  seen  in  Figure  18.  In 
particular,  the  trough  is  placed  immediately  south  of  the  auroral 
oval  and  is  defined  only  from  sunset  to  sunrise,  th,.-  height  of  the 
F2  maximum  ionisation  (h^jF2)  being  50  km  higher  in  the  trough  than 
in  its  absence,  as  indicated  in  Figure  19.  The  location  of  the 
equatorward  boundary  and  the  depth  of  the  trough  can  be  varied  so 
that  fQF2  values  in  the  trough  and  its  latitudinal  extent  can  be 
made  to  agree  with  observations,  when  available.  It  is  noted  that 
Rush  et  al  (Ref.  66)  did  not  substantiate  their  choice  of  50  km  for 
the  trough  height. 


(after  Rush  et  al  (1982) 


Date;  1981/01/20  Time  -  23  ;  00  00 
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In  the  light  of  the  above-mentioned  observations  of  trough  height, 
which  are  summarised  in  Table  3,  it  seems  unfortunate  that  Lockwood 
(Ref.  54)  and  Rush  et  al  (Ref.  65)  chose  a  representative  trough 
height  as  small  as  50  km.  On  the  other  hand,  the  entries  in  Table 
3,  except  for  Lockwood's,  might  well  include  the  more  "spectacular", 
i.e.  easily  identifiable,  examples  of  troughs,  i.e.  include  troughs 
with  greater  than  average  heights. 

The  eibove  discussion  underlines  the  lack  of  information  available 
on  the  height  of  the  mid-latitude  trough.  It  clearly  indicates  the 
desireibility  of  the  continuation  of  observations  of  the  trough  by 
instruments  such  as  the  Advanced  Ionospheric  Sounder  at  Halley  Bay 
(Rodger  and  Pinnock  (Ref.  63);  Dudeney  et  al  (Ref.  23)). 


TABLE  3. 


REFERENCES 

REPRESENTATIVE  MID-LATITUDE 

TROUGH  CHARACTERISTICS 

Height  (Ah) 
km 

Al ti tude 
of  Peak 

km 

Critical 
Frequency  (fQF2) 
MHz 

12.  Buchau  et  al  (1973) 

200 

450 

2 

11.  Bowman  (1969) 

400 

630 

1.5-2. 2 

6  3.  Rodger  and  Pinnock  (1980) 

360,100 

600,430 

2.0,  1.5 

23.  Dudeney  et  al  (1982) 

200 

470 

2 

53,54.  Lockwood  (1980,  1981) 

50 

- 

- 

Nevertheless,  the  information  in  Table  3  is  adequate  to  assign  an 
approximate  value  of  Ah  =  200  km  for  the  trough  height,  which  is  used 
in  the  present  study  in  the  application  of  Halcrow  and  Nisbet's 
(Ref.  39)  model.  Ah  is  the  altitude  of  the  electron  density  rr.ximum 
witliin  the  trough  minus  the  altitude  of  that  value  of  electron  density 
outside  the  trough. 

4.5  Mid-latitude  Trough  Model 

It  is  noted  that  Halcrow  and  Nisbet's  model  provides  contours  of  peak 
F2  electron  density  only.  In  following  ray- tracing  procedures,  it  is 
also  necessary  to  define  electron  density  contours  below  the  peak. 

The  following  guidelines,  additional  to  those  of  Halcrow  and  Nisbet 
(Ref.  39),  were  applied  in  the  present  computations. 

(1)  The  critical  frequency  (f^F2)  within  the  trough  was  obtained  by 
multiplying  the  critical  frequency  in  the  absence  of  the  trough 
by  0.5  (=  if'’).  (it  is  noted  (e.g.  Ref.  9)  that: 

f^  oc  N 

where  f  represents  the  plasma  frequency,  and 

N  the  corresponding  electron  number  density,  so  that: 

[fnF2  (trough  present)*!  ^  _  Nj„F2  (trough  present)  _  ^ 
fQF2 (trough  absent)  J  Nj„F2  (trough  absent)  '  ^ 

The  above  guideline  therefore  corresponds  to  the  Halcrow  and 
Nisbet  procedure  of  multiplying  the  F2  region  electron  number 
density  (NmF2)  in  the  absence  of  the  trough  by  f  =  0.25  to 
obtain  Nj^F2  for  the  trough. 
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(2)  It  is  assumed  that  while  the  trough's  electron  density  is  being 
depleted,  the  altitude  of  its  maximum  is  increasing,  and  becomes 
approximately  200  km  above  the  altitude  at  which  electron  density 
has  the  same  value  outside  the  trough;  i.e.  the  trough  height 
becomes  approximately  Ah  =  200  km.  Due  to  the  reduction  in 
electron  density  within  the  trough,  it  is  seen  from  the  schematic 
diagram  in  figure  20  that  the  actual  altitude  rise  Aa  of  the 
depleted  peak  would  be  less  than  200  km.  For  the  sake  of  compu¬ 
tation  in  Appendix  4,  Aa  was  taken  to  have  a  notional  value  cf 
150  km,  to  sufficient  accuracy  for  the  present  purposes. 

(3)  Figures  9a  and  9b  show  curves  of  constant  plasma  frequency  as 
true  height  versus  latitude  plots  for  specified  simspot  numbers. 
Curves  of  this  kind  represent  the  averaging  of  m2my  troughs  which 
occurred  over  a  wide  range  of  latitudes.  These  troughs  produced 
a  broad  peak  centred  on  the  particular  latitude  (e.g.  57.5  S  in 
Figure  9)  which,  presumably,  was  the  predominant  mean  trough 
position  at  the  corresponding  sunspot  number.  The  height  of 

a  curve,  i.e.  the  altitude  of  a  curve  at  its  peak  {e.g.  at 
latitude  57.5  S)  compared  with  its  altitude  remote  from  that 
latitude,  is,  of  course,  much  less  than  expected  for  individual 
troughs.  However,  the  ratios  of  the  heights  of  the  various 
curves  should  be  the  s«une  as  for  a  representative  individual 
trough.  Absolute  values  for  the  heights  of  curves  that  would 
represent  an  individual  trough  can  be  computed  on  the  assump¬ 
tion,  from  (2)  above,  that  the  height  rise  for  the  curve  repre¬ 
senting  the  critical  frequency  was  such  that  Ah  was  approximately 
200  km. 

In  Appendix  4,  the  above  guidelines,  modified  for  practical  applica¬ 
tion,  are  employed  in  computing  model  mid-latitude  troughs. 

4.6  Planetary  K-index  (Kp)  in  an  Ionospheric  Model 

Figure  21,  derived  from  Ref.  82,  gives  the  probability  distribution, 
over  several  sunspot  cycles,  for  Kp  which  is  a  measure  of  geomagnetic 
activity.  For  instance,  it  indicates  that  the  occurrence  probed^ility 
of  Kp  >  1.5  is  50%,  and  for  Kp  >  3  it  is  81%. 

In  the  present  discussion,  the  occasions  when  the  trough  has  moved 
appreciably  equatorwards  from  the  polar  regions  are  of  most  interest. 

At  such  times,  the  trough  is  more  likely  to  be  located  between  a 
distress  signal  source  in  the  Southern  Ocean  emd  a  receiving  station, 
e.g.  in  Adelaide.  Histograms  denoted  by  asterisks  in  Figure  22 
(adapted  from  Ref.  68)  indicate  the  frequencies  of  occurrence  of  the 
equatorial  boundaries  of  troughs  at  particular  Northern  Hemisphere 
invariant  latitudes.  It  is  assumed  that  similar  invariant  latitudes 
apply  for  the  Southern  Hemisphere,  and  Figure  22  indicates  the  corres¬ 
ponding  southern  geographic  latitudes  (for  the  longitude  of  Adelaide) . 
The  data  in  Figure  22  were  obtained  in  the  interval  February-December , 
1969,  at  a  time  of  maximum  sunspot  activity  when  the  excursions  of 
geomagnetic  activity,  denoted  by  Kp,  would  also  be  maximal. 

Figure  22(b)  indicates  that  even  during  a  year  of  maximum  sunspot 
activity,  there  were  times  of  reduced  geomagnetic  activity  (Kp  <  3+) 
when  the  equatorwards  boundary  of  troughs  would  have  been  over  southern 
Australia,  i.e.  at  geographic  latitudes  equatorwards  of  43  S.  It  may 
be  inferred  from  these  observations  that  during  sunspot  minimum,  when 
Kp  would  normally  have  smaller  values,  there  would  often  be  occasions 
when  the  equatorwards  boundaries  of  troughs  would  occur  at  geographic 
latitudes  of  southern  Australia. 
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Figure  20.  This  schematic  diagram  indicates  that  if  the  peak 
electron  density  (N^)  at  a  given  height  is  reduced 
in  the  formation  of  a  trough,  the  new  peeOc  electron 
density,  e.g.  N|j/4,  would  be  x  km  higher  in  altitude  than 
the  previous  level  at  which  electron  density  was  Nq/4, 
(Outside  the  trough,  electron  density  at  this  (lower) 
level  is  unaltered  at  Nq/4) . 

In  addition,  the  whole  electron  density  profile  is 
raised  ia,  so  that  the  trough  height  becomes  Ah  =  Aa  +  x. 


OauRRENCE  FNEOUENCY  (»WCM)  £  OCCV««na  FNEOUENCV  (pi^) 
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INVARtANT  latitude  (  d*4ra« ) 


Equatorial  Bovmdcury  of  the  Ionization  Irregularity  Zone 
in  the  Northern  Hemisphere  (Kp5.3+)  as  a  Function  of  Local 
Time  for  Low  Magnetic  Activity. 


(b)  Equatorial  Boundary  of  the  Ionization  Irregularity  Zone 

in  the  Northern  Hemisphere  (1^  ^  3+)  as  a  Function  of  Local 
Time  for  High  Magnetic  Activity. 


Figure  22.  Stars  indicate  distributions  which  also  represent 
the  equatorial  boundaries  of  mid-latitude  troughs 
(from  Sagolyn  et  al  (1973),  Ref.  68). 
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It  is  concluded  that  mid-latitude  troughs  can  occur  at  unpredictable 
times  at  sufficiently  low  latitudes  in  the  vicinity  of  Australia  to 
necessitate  making  allowance  for  their  presence  in  conputing  the 
location  in  the  Southern  Ocean  of  the  sources  of  H.F.  ionospherically- 
propagated  radio  signals  received  in  southern  Australia,  e.g.  in 
Adelaide. 

In  enploying  an  ionospheric  model,  which  includes  a  mid-latitude 
F  region  trough,  reliable  values  of  all  parameters  are  readily  avail¬ 
able  except  for  Kp  which  can  vary  rapidly  and  by  large  amounts  in 
response  to  geomagnetic  disturbances  associated  with  solar  activity. 

In  observations  from  September  1962  to  March  1963,  Muldrew  (Ref.  59) 
found  from  73  Alouette  I  topside  soundings  of  troughs  that  unit 
increase  in  Kp  was  associated  with  approximately  3.5  -^eg  equatorward 
movement  of  the  trough  in  geographic  latitude.  (At  the  longitude  of 
Ottawa,  this  corresponds  to  3.4  deg  reduction  in  invariant  latitude.) 

Kohnlein  and  Raitt  (Ref.  46)  enployed  approximately  300  observations 
of  the  trough  obtained  by  means  of  the  ESHO  4  satellite,  from  November 
1972  to  April  1974.  They  determined  an  eitpirical  relationship  for 
the  mean  invariant  latitude  A  of  the  trough  in  terms  of  Kp  and  t 
(local  time  in  hours,  being  positive  after  midnight  and  negative 
before)  namely j 

A  =  65.2  -  2.1Kp  -  0.5t  ±  2  deg 

The  linear  dependence  of  A  on  Kp  is  illustrated  in  Figure  23  (derived 
from  Ref.  46)  which  also  indicates  the  corresponding  south  geographic 
latitudes,  assuming  the  longitude  of  Adelaide. 

Earlier  work  by  Rycroft  and  Burnell  (Ref.  66)  involving  39  topside 
observations  from  Alouette  I,  from  April  to  August  1963,  indicated 
a  relationship: 

A  =  62.7  -  1.4Kp  -  0.7t  ±  3.5  deg 
while  45  observations  by  Rycroft  and  Thomas  (Ref.  67)  indicated: 

A  =  64.6  -  l.SKp  -  0.5t  ±  2.6  deg 

The  above  data  suggest  that  unit  increase  in  Kp  was  associated  with 
decreases  in  invariant  latitude  of  3.4,  2.1,  1.4  and  1.8  deg  respec¬ 
tively,  from  which  2  deg  decrease  is  taken  to  be  representative  of 
these  topside  trough  observations. 

From  October  1964  to  January  1966,  ground  level  observations  of  the 
ionospheric  electron  content  were  made  by  Liszka  (Ref.  52)  employing 
a  Faraday  rotation  technique  at  40  6ind  41  MHz  observing  the  S-66 
satellite.  Some  of  his  observations  are  summarised  in  Figure  24, 
which  indicates  that  a  unit  increase  in  Kp  is  associated  with  a 
reduction  in  geographic  latitude  of  approximately  2.9  deg,  which,  at 
the  longitude  of  Kiruna,  corresponds  to  a  reduction  of  approximately 
3.1  deg  in  invari^ult  latitude. 

On  the  basis  of  much  more  extensive  data,  extending  over  8  years  of 
Alouette  I  amd  II  topside  sounding  observations,  Halcrow  and  Nisbet 
(Ref.  39)  devised  a  trough  model,  examples  of  which  are  shown  in 
Figure  25,  and  which  at  zero  hours  local  time  is  represented  as 
follows: 


GEOINAGNETIC  ACTIVITY 


41 


ERL-0449-SD 


South  GoograpMc  Lalitudo  at 


(a)  Position  of  the  Mid-latitude 
Trough:  Invariant  Latitude  vs 
Geomagnetic  Activity  at  19.6  hr  LT. 

The  dots  refer  to  the  Northern  Hemis¬ 
phere  and  the  circles  to  the  Southern 
Hemisphere. 


(b)  Position  of  the  Mid-latitude 
Trough:  Invariant  Latitude  vs 
Geomagnetic  Activity  at  1 . 2  hr  LT 
(Northern  Hemisphere) . 


(c)  Decrease  of  Invariant  Latitude 
with  Increasing  Geomagnetic 
Activity. 

t  Gradient  at  different  local 
times. 

— Average  gradient  of  all  local 
times  investigated. 

KOHNLEIN  AND  RAITT  (1977) ,  Ref.  46 


Figure  23. 
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Figtire  25(c)  .  Modified  CCIR  Predicted 
Nj^Fj  in  m"^  for  Spring 
Equinox;  =  50,  Kp  =  0+ 
(after  Halcrow  and  Nisbet 
(1977) .  Ref.  39) . 


Figure  25  (d)  .  Modified  CCIR  Predicted  in  m" 

Summer  Solstice;  =  50,  Kp  =  0+ 
(after  Halcrow  and  Nisbet  (1977), 
Ref.  39) . 


Winter  Solstice;  R^  =  50,  Kp  =  0+ 
(after  Halcrow  and  Nisbet  (1977), 
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0.5(K.p  -  i) 
1.0(ltp  -  j) 
1.25(Kp  -  i) 


where  A  represents  invciriant  latitude,  the  subscripts  T,  B,  S  and 
N  represent  top,  bottom,  equatorward  and  poleward  respectively, 
and  the  non-fractional  numbers  represent  degrees. 


In  particular,  these  formulae  indicate  a  decrease  of  the  order  of 
1  deg  in  invariant  latitude  per  unit  increase  in  Kp,  which  is  a 
smaller  decrease  than  those  indicated  from  the  earlier  topside  and 
bottomside  observations  mentioned  above. 


In  the  light  of  the  above  observations,  it  appears  that  Halcrow 
and  Nisbet's  (Ref.  39)  model  may  well  be  improvable  as  regards  the 
magnitude  of  the  dependence  of  A  on  Kp.  It  highlights  the  impor¬ 
tance  of  mcJcing  direct  observations  of  the  position  of  the  mid- lati¬ 
tude  trough,  prefer^d5ly  in  as  nezir  to  real  time  as  possible  for 
emergency  procedures,  rather  than  relying  too  heavily  on  a  model 
that,  of  necessity,  incorporates  average  values  of  the  trough  posi¬ 
tion.  It  also  underlines  the  desiredsility  of  obtaining  statistical 
data  for  the  trough  position  in  the  Southern  Hemisphere. 

The  question  of  how  rapidly  the  ionosphere  responds  to  a  change  in 
Kp  requires  consideration.  Topside  sounder  data  presented  in  Refs. 

59  and  36  suggested  that  an  increase  in  Kp  preceded  a  reduction  in 
the  trough  latitude  by  several  hours,  namely  of  the  order  of  5  hours 
and  17  hours  respectively.  These  rough  estimates  are  in  qualitative 
agreement  with  Feinblum  and  Horan's  (Ref.  32)  more  general  conclusion 
that  the  F2  region  critical  frequency  (foF2)  was  best  correlated 
with  Kp  observed  some  12  hours  earlier. 

In  the  present  context,  the  several  hours'  delay  between  Kp  and  the 
F2  region  response  is  of  academic  interest  only.  Information  from 
the  Geophysikalisches  Institut  at  the  University  of  Gottingen, 

Federal  Republic  of  Germany,  which  provides  an  International  Service 
of  Geomagnetic  Indices  based  on  observations  at  many  locations , 
indicated  that  there  are  4  to  6  weeks'  delay  between  geomagnetic 
observations  being  made  and  the  derived  Kp  index  being  available. 

(Kp  computations  are  based  on  local  K  (i.e.  magnetic  activity)  obser¬ 
vations  at  13  stations,  including  Canberra  in  Australia  and  Eyrewell 
in  New  Zealand.) 


The  non-availability  of  Kp  values  in  something  approaching  real  time 
presented  a  dilemma.  A  possible  solution  was  to  obtain  K  values  of 
geomagnetic  variations  from  Australian  observing  stations,  in  order 
to  compute  a  parameter  compareible  with  Kp.  Problems  were  anticipated 
in  co-ordinating  such  a  procedure,  which  would  involve  attempting, 
often  outside  of  normal  working  hours,  to  obtain  data  by  long-distance 
communication  with  observers.  This  suggested  that  a  preferable  alter¬ 
native  would  be  to  observe  the  position  of  the  mid-latitude  trough 
directly. 

A  direct  and  convenient  technique  for  locating  the  position  of  the 
trough  would  be  to  record,  e.g.  adjacent  to  a  direction-of-arrival 
measuring  station  in  Adelaide,  the  Doppler  shift  of  dual-frequency 
H.F.  transmissions  from  a  circumpolar  orbiting  satellite.  In  this 
way,  changes  in  total  electron  content  along  the  continually  changing 
satellite-to-receiver  path  could  be  determined,  and  the  trough  position 
identified.  It  would  be  desirable  to  use  several  such  satellites, 
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so  that  the  trough  position  could  be  monitored  at  intervals  of  no 
more  than  an  hour  or  so.  It  would  also  be  desirable  for  these 
satellites  to  transmit  on  the  same  pair  of  coherent  frequencies. 

There  are  five  U.S.  Navy  Navigational  Satellites  which  satisfy 
this  criterion  and  transmit  coherently  at  150  and  400  MHz  in  circum¬ 
polar  orbits.  The  dual-frequency  Doppler  technique  of  measuring 
ionospheric  electron  content  is  described  and  discussed  in  Section  5. 

Another  requirement  of  a  dual-frequency  Doppler  receiving  facility 
is  that  it  be  immediately  available  for  determining  the  trough 
location  at  times  of  emergency.  Unfortunately,  the  existing  TRANET 
satellite  tracking  facility  at  Smithfield,  near  Adelaide,  does  not 
satisfy  this  latter  requirement,  since  its  use  is  dedicated  to 
United  States  Government  purposes.  However,  dual  frequency  TRANET 
data  would  be  available  for  statistical  studies  of  the  mid-latitude 
trough  position,  through  courtesy  of  the  U.S.  Defense  Mapping  Agency. 

4.7  Scintillation  Fading 

Lyon  (Ref.  55)  has  indicated  the  potentially  adverse  effects  of 
auroral  ionospheric  scintillations  on  Doppler  frequency  measurements 
using  NNSS  satellites,  discussed  in  Section  5.3  below.  In  particular, 
to  ensure  continuous  signal  reception,  a  good  phase-locked  receiving 
system  is  required  which  will  hold  the  signal  through  severe  scintil¬ 
lations,  allowing  the  fluctuations  to  be  smoothed  with  the  use  of 
suitable  time  constants. 

Similar  adverse  effects  are  expected  with  the  proposed  SARSAT  Satellite 
system  (discussed  in  Section  5.2  below)  which  will  receive  VHF 
(121.5  MHz)  distress  transmissions  through  the  ionosphere.  At  higher 
geomagnetic  latitudes,  scintillations  due  to  ionospheric  irregulari¬ 
ties  will  degrade  the  signal  by  imposing  rapid  time  fluctuations, 
an  effect  that  is  reduced  at  higher  operating  frequencies.  It  is 
concluded  in  Ref.  55  that  in  moderate  scintillating  conditions,  a 
VHF  signal  will  be  acquired  without  large  errors  in  estimating  the 
position  of  a  distress  signal  source;  but  in  severe  scintillation 
conditions,  the  signal  will  simply  not  be  acquired  at  high  geomag¬ 
netic  latitudes.  Table  4  includes  data  from  Ref.  14  which  emphasises 
the  severity  of  scintillation  fading  at  auroral  latitudes,  particu¬ 
larly  at  night  under  magnetically  disturbed  (high  Kp)  conditions. 

The  present  project  would  afford  an  opportunity  to  gain  additional 
data  on  scintillations  at  high  invariant  geomagnetic  latitudes,  as 
measured  by  the  fading  of  NNSS  satellite  transmissions  at  150  and 
400  MHz. 

4.8  Ground-Based  Transmitting  Stations 

Further  scintillation  data  are  obtainable  by  recording  at  Adelaide, 
fixed  frequency  HF  transmissions,  e.g.  from  Mawson  and  Casey  (by 
courtesy  of  the  Australian  Ionospheric  Prediction  Service)  and 
Macquarie  Island.  However,  a  more  important  purpose  in  receiving 
such  transmissions  is  to  chec)t  their  observed  directions  of  arrival 
against  the  computed  values.  In  short,  a  valuedsle  check  is  provided 
on  the  suitability  of  Halcrow  and  Nisbet's  (Ref.  39)  Northern  Hemis¬ 
phere  raid-latitude  trough  model  as  an  addition  to  the  RADFIN  program 
of  DRCS  for  Southern  Hemisphere  studies. 

The  present  project  could  involve  only  one  direction  of  arrival 
aerial  system,  which,  in  principle,  should  be  sufficient.  However, 
experiments  in  the  reception  of  HF  transmissions  from  ground-based 
stations  may  indicate  that  two  or  more  such  systems  in  southern 
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Australia  are  desirable  to  ensure  (i)  sufficient  accuracy  in  locating 
the  position  of  a  distress  signal  source,  and  (ii)  that  adeauate 
signal  strengths  are  likely  to  be  received  by  at  least  one  system 
under  conditions  of  severe  ionospheric  disturbance. 

TABLE  4.  PERCENTAGE  OF  OCCURRENCE  OF  SCINTILLATIONS* 


(a)  12  dB  pe«Uc  to  peak  at  137  MHz,  sub-auroral  and 

auroral  latitudes 


SITE 

Kp 

DAY 

(0500*1700  LT) 

NIGHT 

(1700-0500  LT) 

Sagamore  Hill  (Massachusetts) 

0  to  3+ 

0 

1.4 

>  3+ 

0.1 

2 

Goose  Bay  (L«ibrador) 

0  to  3+ 

>  3+ 

0.1 

1.6 

1.8 

6.8 

Narssarssuaq  (Greenland) 

0  to  3+ 

18 

45 

(b)  ^  10  dB  pecik  to  peak  at  254  MHz,  auroral  latitudes 


SITE 

DAY 

NIGHT 

(0600-1800  LT) 

(1800-0600  LT) 

Goose  Bay  (Labrador) 

0 

to  3+ 

0.1 

0.1 

> 

3+ 

0.3 

1.2 

Narssarssuaq  (Greenland) 

0 

0.1 

0.9 

> 

3+ 

2.6 

8.4 

LT:  Local  Time  *  (  ata  from  Ref.  14) 


5.  SATELLITES  -  SEARCH  AND  RESCUE  CAPABILITIES 
5.1  Geostationary  Satellites 

There  are  currently  available  several  communications  satellites 
in  geostationary  orbits  (over  the  equator)  with  the  capacity  to 
provide  assistance  in  recording  distress  signals,  e.g.  three 
geostationary  MARISAT  satellites  (Ref.  20)  are  located  over  the 
Atlemtic  Ocean  (15  W) ,  the  Indian  Ocean  (a  "spare"  satellite,  73  E) 
and  the  Western  Pacific  (176.5  E) .  These  provide  coiranercial  communi¬ 
cations  services,  via  satellite,  between  ships  and  shore  stations 
(at  Southbury,  Connecticut,  for  the  Atlantic  satellite,  and  Santa 
Paula,  California,  for  the  Pacific  satellite)  and  thence  to  existing 
public  teleprinter,  telephone  and  facsimile  communications  systems. 
Figure  26  shows  the  service  coverage  of  the  MARISAT  system,  assuming 
elevation  angles  exceeding  5  deg  from  ship  to  satellite.  The  coverage 
of  the  Indian  Ocecin  and  Pacific  satellites  is  seen  to  include  most 
of  the  Southern  Ocean  south  and  south-west  of  Australia,  which  is 
of  present  concern. 

The  configuration  of  the  MARISAT  Communications  System  incorporates, 
as  an  objective,  the  provision  of  "improved  safety  and  distress 
services",  ouid  permits,  as  a  feature,  the  "handling  of  emergency 
calls  on  a  priority  basis  over  other  traffic".  However,  the  expense 
of  the  comparatively  sophisticated  ship-board  equipment  on  L-band , 
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at  1638.5-1642.5  MHz,  for  transmission,  and  1537-1541  MHz,  for 
reception,  together  with  a  steerable  antenna  to  be  directed  towards 
the  satellite,  would  preclude  such  installations  from  any  but  the 
best-equipped  ships,  e.g.  many  overseas,  and  some  Australian, 
tankers. 

Furthermore,  it  is  expected  that  ships  using  the  MARISAT  system 
will  provide  regular  updates  of  their  exact  positions  to  shore 
stations.  The  location  of  a  ship  in  distress  could  therefore  be 
accurately  estimated  from  a  previous  position  without  the  need  to 
take  bearings  on  it.  Transmit/receive  aerials  on  the  MARISAT 
satellites  are  "despun  and  earth  oriented",  i.e.  stationary  relative 
to  the  earth's  surface;  thus,  even  though  the  satellites  could  detect 
an  "unknown"  distress  signal  on  the  appropriate  frequency,  its  bearing 
could  not  be  determined. 

Each  MARISAT  satellite  also  has  three  ultra-high  frequency  (UHF) 
channels  in  the  250-300  MHz  band  reserved  for  U.S.  Navy  use;  these 
would  be  subject  to  some  adverse  ionospheric  effects  in  the  form 
of  scintillation  fading.  The  fact  that  these  frequencies  are  similar 
to  the  VHF  international  distress  frequencies,  121.5,  156.8  and 
243  MHz  indicates  the  suitability,  in  principle,  of  geostationary 
satellites  such  as  MARISAT  to  record  VHF  distress  signals. 

Ship  and  'plane  communication  has  traditionally  been  below  250  MHz, 
i.e.  at  VHF  and  lower.  'Planes,  and  the  majority  of  shipping,  would 
not  be  equipped  to  exploit,  in  times  of  emergency,  GHz  frequencies 
available  commercially  with  MARISAT.  Nevertheless,  the  advantages 
of  the  wide  coverage  of  the  earth's  surface  and  the  lack  of  gross 
ionospheric  fading  effects  at  GHz  frequencies  supports  the  trend 
to  their  increasing  usage  in  geostationary  satellites. 

The  MARISAT  system  was  described  above  as  an  example  of  a  geosta¬ 
tionary  satellite  system;  its  advantaaes  and  limitations  in  land/sea 
rescue  procedures  were  indicated. 

The  INMARSAT  system,  which  is  intended  to  take  the  place  of  MARISAT, 
will  involve  the  use  of  six  satellites  (orbiting  and  stationary) 
and  include  the  facility  for  real  time  sea  distress  radio  transmis¬ 
sions  (Ref.  26) . 

5.2  Orbiting  Satellites  -  Proposed  SARSAT  System 

It  is  reported  (Ref.  26)  that  the  U.S. A.,  Canada  and  France  are 
currently  co-operating  on  a  proposed  SARSAT  (Search  and  Rescue 
Satellite  System)  which  will  employ  three  weather  satellites  of 
the  TIROS-N  (Television  Infrared  Radiation  Observation  Satellite) 
system.  These  will  have  almost  circular  near-polar  orbits  with  a 
98.7  deg  inclination  to  the  equatorial  plane  at  a  height  of  834  km, 
giving  a  period  of  101.5  min.  The  satellites  will  each  be  equipped 
with  a  SARSAT  transponder,  and  utilise,  initially,  emergency  position 
indicating  radio  beacons  (EPIRB)  operating  at  the  121.5  and  243  MHz 
international  distress  frequencies  at  transmitted  powers  of  the 
order  of  100  mW.  With  this  system,  positions  of  EPIRB' s  are  claimed 
to  be  determinable  within  10  to  20  km.  It  is  proposed  that,  later, 
beacons  (being  developed)  will  be  employed  at  406  MHz  and  allow  posi¬ 
tions  to  be  determined  within  2  to  5  km.  In  all  cases  the  emergency 
transmission  is  received  at  the  satellite  in  an  on-board  repeater, 
translated  to  a  kHz  band  of  frequencies  which  is  transmitted  to  a 
ground  station  as  phase  modulation  on  a  1544.5  MHz  down-link.  Because 
of  the  satellite  motion,  the  received  emergency  transmissions  are 
Doppler  shifted,  and  the  information  on  emergency  transmitter  location 
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is  contained  in  the  Doppler  signature  in  the  usual  manner  for  a 
Doppler  navigation  system,  such  as  the  Transit  or  Global  Positioning 
Systems  mentioned  in  Sections  5.3  and  5.4  below.  Data  processing 
for  position  location  is  carried  out  at  the  ground  stations. 

(Appendix  2  summarises  the  techniques  of  Doppler  measurement  of  range.) 

Lyon  (Ref.  55)  has  shown  that  in  locating  the  emergency  transmitter 
at  a  single  frequency  of  121.5  MHz,  errors  as  great  as  35  km  in 
ground  range  can  arise  due  to  horizontal  gradients  in  the  mid-latitude 
trough  (described  in  Section  6)  or  high  latitude  trough  (Section  6.1) 
in  F  region  electron  density.  The  ionospheric  contribution  to  the 
Doppler  frequency  change  is,  however,  inversely  proportional  to 
frequency  so  that  errors  arising  from  the  Doppler  effect  would  be 
correspondingly  reduced  at  a  higher  frequency,  e.g.  406  MHz.  Lyon 
(Ref.  55)  also  indicated  that  at  high  latitudes,  rapid  fading 
(scintillations)  due  to  ionospheric  irregularities  will  degrade 
the  emergency  signal.  He  suggested  that  with  severe  scintillations, 
the  signal  will  not  be  acquired  at  high  latitudes.  (Clearly,  EPIRB 
transmissions  at  two  frequencies,  e.g.  121.5  and  406  MHz  could  be 
employed,  in  principle,  to  locate  EPIRB' s  to  appreciably  greater 
accuracy  by  employing  a  differential  Doppler  technique  analogous  to 
that  described  in  Chapter  5.) 

5.3  Orbiting  Satellites  -  NNSS 

The  Navy  Navigational  Satellite  System  (NNSS)  or  Transit  system, 
also  called  the  NAVSAT  system,  has  been  operational  since  1964, 
available  for  commercial  use  since  1967  and  is  expected  to  provide 
valuable  service  as  a  navigation  aid  for  shipping  until  at  least 
1995,  after  which  phase-over  to  a  Global  Positioning  System  (GPS) 
of  navigation  (described  in  Section  5.4  below)  is  expected  to  be 
complete  (Ref.  79) .  The  NNSS  is  operated  by  the  Navy  Astronautics 
Group  with  headquarters  at  Point  Magu,  California,  and  is  under 
subsidy  from  the  United  States  Navy.  At  any  time  there  are  a  maxi¬ 
mum  of  6  (currently  5)  operational  satellites  in  circular  polar 
orbits  at  about  1075  km  altitude  and  having  periods  of  107  min. 

The  orbits  form  a  "bird-cage"  (Figure  27)  within  which  the  earth 
rotates  carrying  an  earth  observer  past  each  orbit  in  turn.  In 
principle,  an  observer  can  obtain  a  position  "fix"  whenever  a 
satellite  is  above  the  horizon.  (The  Doppler  technique  for  position¬ 
fixing  is  outlined  in  Appendix  2.) 

In  practice,  passes  are  accepted  between  8  and  70  deg  of  elevation. 
Tropospheric  refraction  produces  excessive  errors  (over-estimates) 
in  range  below  8  deg  (Stansell  (Ref.  76):  his  Figure  44).  At  eleva¬ 
tions  cibove  70  deg,  limitations  in  knowledge  of  the  observer's 
altitude  above  the  geoid  can  result  in  excessive  error  in  a  "fix" 
of  the  observer's  longitude  (Stansell  (Ref.  76);  his  Figure  46). 

This  longitude  error  arises  as  follows  for  a  satellite  in  polar 
orbit.  Figure  28  depicts  a  satellite  moving  perpendicular  to  the 
page  and  at  its  closest  approach  to  stations  X,  Y  and  Z  at  alti¬ 
tude  1.  Stations  X',  Y'  and  Z'  at  altitude  2  would  have  respectively 
the  same  rcinges  and  rates  of  change  of  range  as  for  X,  Y  and  Z.  In 
particular,  both  longitude  amd  altitude  affect  the  range  and  its 
rate  of  change  in  similar  ways.  In  computing  an  accurate  longitude 
fix,  it  is  therefore  necessary  to  have  prior  knowledge  of  the 
receiving  station's  altitude. 

With  five  satellites,  the  average  time  interval  between  fixes  varies 
from  about  32  min  at  70  deg  latitude  to  about  97  min  at  the  equator, 
and  is  81  min  at  35  deg  latitude  (Stansell  (Ref.  77);  his  Figure  32). 
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Fiqure  28.  Effect  of  Altitude  Estimate  on  Position  Fix  (after 
Stansell  (1978),  Ref.  77). 
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In  the  context  of  the  present  report,  NNSS  satellites  would  be 
used  for  finding  the  real  time  position  of  the  F  region  raid- latitude 
trough.  This  would  involve  differential  Doppler  measureraents  of 
ionospheric  electron  density  using  150  and  400  MHz  transmissions 
at  powers  of  approximately  J  W  and  1  W  respectively.  Prior  know¬ 
ledge  is  assumed  of  the  location  of  a  receiving  station,  e.g.  near 
Adelaide,  and  of  the  particular  received  satellite  at  any  instant. 

5.4  Orbiting  Satellites  -  Global  Positioning  System 

Currently,  the  Global  Positioning  System  (GPS)  of  navigation  involves 
six  satellites. 

The  GPS  NAVSTAR  system  will  be  a  new  satellite  constellation 
consisting  of  18  satellites  distributed  uniformly  within  six  orbital 
planes,  each  at  55  deg  inclination  to  the  equator;  the  first  satellite 
is  scheduled  for  launching  in  1983,  the  eighteenth  in  1988.  The 
satellites  orbit  at  about  17.400  km  above  the  earth,  completing  an 
orbit  cycle  every  12  hours.  They  use  frequencies  of  1575.42  and 
1227.6  MHz  (Refs.  76  and  10). 

The  objective  of  NAVSTAR  is  to  provide  direct  line-of-sight  signals 
continuously  from  at  least  four  satellites  to  all  users  on  a  world¬ 
wide  basis.  A  ground  control  network  tracks  the  satellite  signals, 
determines  the  orbits  precisely,  and  transmits  an  appropriate  orbit 
definition  to  the  memory  system  of  each  satellite.  By  this  means,  the 
satellites  can  provide  users  with  an  accurate  description  of  their 
positions  as  a  function  of  time.  To  obtain  a  position  fix,  the  user 
must  be  equipped  with  a  receiver  capable  of  tracking  signals  from 
four  satellites,  either  simultaneously  or  sequentially.  Each  set 
of  four  time-of-arr:.val  measurements  of  signals  made  by  the  user 
permits  four  independent  variables  to  be  determined,  namely  clock 
time  offset,  latitude,  longitude  and  altitude.  Hence,  the  following 
accuracies  are  achievable:  position,  16  m;  velocity,  0.1  ms”'; 
time,  100  ns.  GPS  provides  near  real-time  determinations  of  three- 
dimensional  position  and  velocity:  (the  TRANSIT  system,  described 
in  Section  5.3  above,  permits  two-dimensional  determinations  only). 

A  future  phase  proposed  for  the  GPS  system  involves  an  additional 
28  satellites  in  three  orbital  planes.  These  would  be  Space  Shuttle 
launched,  commencing  in  1986. 

In  the  present  context,  the  frequencies  (1575.42  and  1227.6  MHz) 
employed  by  GPS  would  be  too  high  to  provide  differential  Doppler 
measurements  of  electron  density  with  accuracy  comparable  to  that 
achievable  using  the  TRANSIT  system.  Nevertheless,  the  55  deg 
inclination  of  NAVSTAR  satellite  orbits  means  that  they  would  pass 
over  geographic  latitudes  as  high  as  55  deg  which,  at  the  longitude 
of  Adelaide,  corresponds  to  68  deg  south  invariant  latitude:  see 
Figure  3.  The  GPS  satellites  therefore  traverse  latitudes  suffici- 
ciently  high  that  satellite  transmissions  received  in  Australia 
could  pass  through  the  mid-latitude  trough  under  moderately  disturbed 
geomagnetic  conditions. 

The  GPS  system  would  therefore  provide  the  potential  for  real-time 
measurements  of  the  total  electron  density  in  the  trough,  but  with 
reduced  accuracy  due  to  the  higher  (GHz)  frequencies  employed. 
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6.  IONOSPHERIC  F  REGION  MID-LATITUDE  TROUGH 

6 . 1  General  Features 

The  mid-latitude  electron  density  trough  was  discovered  by  Muldrew 
(Ref.  59),  and  is  a  major  but  incompletely  understood  feature  of 
the  sub-auroral  ionosphere.  It  has  very  significant  effects  on 
radio-wave  propagation  (Refs,  81,  39  and  23) .  It  is  most  prominent 
in  the  winter  night-time  ionosphere,  and  consists  of  a  circumpolar 
depletion  of  the  plasma  concentration  (tlie  "trough  minimum")  corres¬ 
ponding  to  electron  densities  at  the  F2  region  maximum  which  are 
less  by  an  order  of  magnitude  than  maxima  elsewhere. 

Although  the  present  concern  is  m2u.nly  with  the  mid-latitude  trough, 
it  should  be  noted  that  there  is  also  a  polar  trough  v,hich  is  a 
feature  of  very  high  invariant  latitudes  (defined  in  Appendix  1) 
greater  than  about  80  deg  (e.g.  Ref.  23).  Both  troughs  are  most 
clearly  evident  on  winter  nights.  At  other  seasons,  solar  photo¬ 
ionisation  tends  to  reduce  their  prominence. 

6.2  Observations  of  Location  and  Extent  of  the  Trough 

(i)  Topside  Soundings  and  other  Satellite  Techniques 

Since  Muldrew's  (Ref,  59)  initial  description  of  the  "ionos¬ 
pheric  trough"  appeared,  based  on  topside  soundings  from  the 
Alouette  I  satellite,  it  has  been  extensively  studied  by 
satellite  techniques  in  attempting  to  understand  its  spatial 
and  temporal  variations,  e.g.  Refs.  18,  67  and  79.  The  main 
data  gathering  instruments  were  satellite-borne  ion  traps 
(e.g.  Ref.  73)  and  electrostatic  probes  (e.g.  Refs.  80  and 
68)  to  measure  the  ion  density,  ion  mass  and  the  ratio  of 
ion  mass  to  ion  temperature,  as  well  as  topside  ionosonde 
sounding.  Sagolyn  et  al's  (Ref.  6C)  observations  indicate 
that  mid-latitude  troughs  occurred  on  27%  of  all  orbits  at 
all  times;  furthermore,  their  equatorward  edges  had  mean 
invariant  latitudes  of  some  56  to  72  deg  which  at  Adelaide's 
longitude  (138.6  E)  would  correspond  to  mean  geographic  lati¬ 
tudes  of  44  to  54  S.  They  also  showed  that  either  troughs 
or  single  gradient  depressions  (i.e.  appearing  at  the 
equatorial  boundary  but  not  followed  by  a  poleward  increase) 
occurred  on  51%  of  all  orbits  at  all  times.  Figure  29  shows 
a  significant  diurnal  variation  in  the  frequency  of  trough 
occurrence.  Troughs  were  rarely  observed  in  daylight  between 
6  and  16  hrs  local  time;  maximum  occurrence  was  around  4  hrs 
local  time. 

From  observations  of  approximately  3500  orbits  of  Isis  1, 
over  some  12  months,  Sagolyn  et  al  (Ref.  68)  found  that  the 
scale  sizes  of  mid-latitude  troughs  typically  ranged  from 
140  to  1500  km  in  latitudinal  extent;  some  examples  are  shown 
in  the  top  diagrams  of  Figure  29.  The  latitudinal  extent  of 
mid-latitude  troughs  is  also  indicated  in  the  models  in 
Figures  10  and  24  (Ref.  39)  and  Figure  30  (Ref.  54)  discussed 
previously  in  Section  2.4.  The  trough  can  extend  through 
almost  as  much  as  180  deg  of  longitude.  The  overall  invariant 
latitude  width  of  the  reduced  electron  density  region  can  be 
about  10  deg  including  the  trough  edaes. 
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Figure  29  (a)  .  Exeutioles  of  loniaatAon 
Troughs  (Sagolyn  et  al 
(1973),  Ref.  68). 
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Figure  29(b).  A  Transpolar  Orbit 

Showing  Nightside  and 
Dayside  Trough  (Sagolyn 
et  al  (1973) ,  Ref.  68) . 
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(Sagolyn  et  al  (1973), 
Re'.  68). 
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Figure  30 (a) .  Invariant  Latitude  A,  Local  Time  (LT)  Plot  of  the 
Positions  of  the  Predicted  Trough  Walls  for  the 
Night  of  May  18-19,  1976  (Kp  =  1.1). 

(Loc)twood  (1981),  Ref.  54) 


Figure  30(b).  A-LT  Plot  of  Mean  Trough  for  Kp  =  2.75,  May  1976. 
(Loc)wood  (1981),  Ref.  54) 
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The  mid-latitude  trough  is  bounded  on  the  high  latitude  side 
by  a  steep  "wall"  of  ionisation,  the  poleward  "edge",  which 
is  approximately  2-4  deg  of  latitude  wide,  and  is  created 
primarily  by  soft  {<1  keV)  electron  precipitation  from  the 
magnetosphere  (Ref.  62).  in  the  low  latitude  boundary  (the 
"equatorial  edge")  which  is  approximately  3-8  deg  of  latitude 
wide,  the  electron  concentration  is  highly  variable  and  is 
an  extension  of  the  mid-latitude  ionosphere  (Ref.  24).  The 
mean  position  of  the  trough  is  at  about  60  deg  of  invariant 
latitude  (which  gives  a  good  approximation  to  the  actual 
geomagnetic  field,  as  discussed  in  Appendix  1) .  The  trough 
is  most  pronounced  at  the  maximum  of  the  F2  layer  ionisation, 
but  has  been  detected  up  to  2500  km. 

In  the  particular  case  of  night-time  in  winter,  when  the 
trough  was  most  pronounced,  Mendillo  and  ChtOco  (Ref.  57) 
studied  Northern  Hemisphere  topside  ionospheric  sounder 
observations  obtained  with  the  Isis  2  satellite  during  a 
period  of  very  low  magnetic  activity.  They  showed  (Figure 
31)  that  at  the  altitude  of  maximum  electron  density,  the 
poleward  wall  of  the  trough  could  be  represented  as  having 
an  invariant  latitudinal  width  of  about  2  deg,  a  trough 
minimum  approximately  13  deg  wide  and  an  equatorward  edge 
of  the  order  of  7  deg  width,  i.e.  a  structure  (at  a  geographic 
longitude  of  some  70  W)  having  an  overall  invariant  latitude 
width  of  the  order  of  22  deg.  (Mendillo  cuid  Chako's  (Ref.  57) 
data  were  given  in  corrected  geomagnetic  latitude,  considered 
briefly  in  Appendix  1) . 

Feinblum  and  Horan  (Ref.  32)  observed  the  trough  width  to  be 
a  variable  feature  depending  on  local  time,  the  trough  being 
narrow  near  midnight,  and  wide  near  dawn  and  dusk.  Further¬ 
more,  the  trough  width  decreased  as  the  planetary  magnetic 
activity  index,  Kp,  increased.  The  width  of  the  equatorward 
wall,  in  particular,  decreased  with  increasing  Kp.  The 
depth  of  the  trough  also  varied  with  local  time  and  Kp; 
typically,  electron  densities  in  the  trough  were  a  factor 
of  about  5-10  lower  than  those  found  on  either  the  poleward 
or  equatorward  side  of  the  trough.  The  trough  tended  to  be 
deepest  near  midnight,  with  the  depression  in  electron  density 
being  about  a  factor  of  8  to  10.  Past  midnight,  the  trough 
filled  in  slightly,  cund  the  walls  became  less  steep.  Near 
dawn  and  dusk  the  trough  became  less  pronounced  owing  to 
photoionisation . 

(ii)  Bottoms ide  Soundings 

Bottomside  ionogram  observations  of  the  trough  are  rarely 
reported.  This  is  because  a  trough  that  is  directly  eibove 
a  transmitter  will  return  off-vertical  reflections  from  its 
steep  "sides"  (contours  of  equal  ionisation)  in  directions 
some  tens  of  degrees  off-vertical  (Figure  2)  with  critical 
frequencies  close  to  those  that  would  prevail  in  the  absence 
of  the  trough  (Bowman  (Ref.  11)).  As  these  reflections 
produce  ionograms  with  critical  frequencies  and  ranges 
(200-300  km)  similar  to  those  prevailing  in  the  absence  of 
the  trough,  they  are  readily  mistaken  for  normal  F2  region 
echoes.  However,  since  the  troughs  are  deep,  as  seen  in 
Figure  2,  the  vertically  incident  signals  produce  traces 
with  ranges  around  500  km  and  critical  frequencies  several 
times  smaller  than  those  of  the  lower  traces;  vertical 
reflections  from  the  trough  therefore  tend  to  be  "lost" 
amongst  other  reflections.  As  a  result  of  such  effects. 
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only  an  exceptionally  wide  trough  cam  be  resolved  by  means 
of  usual  ionosonde  techniques.  <This  corresponds  to  very 
quiet  geomagnetic  conditions,  or  else  a  location  where  the 
invariant  latitude  lines  are  exceptionally  far  apart . ) 
Alternatively,  special  techniques  must  be  used  to  identify 
the  bottomside  trough;  e.g.  Bowman  (Ref.  11)  employed  an 
ionosonde  with  direction-finding  capabilities. 

Rodger  and  Pinnock  (Ref.  63)  have  identified  five  types  of 
trough  on  ionograms  recorded  with  a  normal  ionosonde  at 
Halley  Bay.  Their  classifications  are  associated  with 
time  of  occurrence,  Kp  and  relationship,  or  otherwise,  with 
sporadic  E  and/or  spread  F  ionisation.  They  made  the  inter¬ 
esting  observation  that  the  level  of  magnetic  activity, 
measured  by  Kp,  required  for  a  particular  trough  sequence 
to  be  observed  at  Halley  Bay  was  slightly  lower  at  sunspot 
maximum  compared  with  sunspot  minimum.  They  suggested  that, 
in  principle,  if  Kp  is  known,  the  time  of  occurrence  of  the 
trough  (e.g.  at  Halley  Bay),  its  shape  and  movement  can  be 
predicted.  This  could  lead  to  a  significant  improvement  in 
the  prediotedjility  of  radio-wave  propagation  in  the  vicinity 
of  the  trough. 

Preliminary  observations  at  Halley  Bay  with  an  Advanced  Iono¬ 
spheric  Sounder  have  demonstrated  its  suitability  for  recor¬ 
ding  the  trough  and  tracking  its  diurnal  motion  by  locating 
its  echoing  region  in  range,  azimuth  and  elevation  (Dudeney 
et  al  (Ref.  23)). 

The  mid-latitude  trough,  along  with  the  auroral  oval,  is 
asymmetrically  located  with  respect  to  the  geographic  polar 
axis,  as  suggested  in  Figure  32,  but  is  fixed  in  space  rela¬ 
tive  to  the  sun.  Consequently,  as  the  earth  rotates,  a 
fixed  observer  on  the  earth's  surface  records  a  diurnal 
motion  of  the  trough  which  moves  equatorwards  in  the  early 
evening  hours  (local  time),  and  afterwards  moves  polewards. 
Trough  velocities  relative  to  a  fixed  observer  2u:e  of  the 
order  of  50-150  m/s  (Bowman  (Ref.  11)).  The  mid-latitude 
trough  therefore  assumes  the  appearance  of  a  very  large- 
scale  travelling  ionospheric  disturbance. 

Lietinger  et  al  (Ref.  51)  made  differential  Doppler  measure¬ 
ments  of  transmissions  from  Navigation  satellites  at  four 
middle  to  high  latitude  European  stations.  During  some 
winter  nights  they  observed  F  region  ionospheric  troughs 
(indicated  in  Figure  33)  with  poleward  boundaries  near  the 
equatorward  edge  of  diffuse  aurorae  coinciding  with  the 
precipitation  of  soft  electrons  (<400  eV) .  They  commented 
that  these  electron  content  observations  permitted  F  layer 
investigations  to  be  made  under  disturbed  conditions  when 
ionograms  failed  or  were  very  difficult  to  interpret. 

In  making  differential  Doppler  measurements  of  NNSS  satellite 
transmissions  at  150  and  400  MHz,  Evans  and  Wand  (Ref.  31) 
observed  that  the  trough  did  occur  in  the  daytime,  although 
it  was  mainly  a  night-time  phenomenon.  They  found  that  on 
some  magnetically  disturbed  winter  days ,  deep  troughs  occurred 
in  the  afternoon  and  evening.  These  troughs  are  exemplified 
in  Figure  34,  and  gave  rise  to  the  largest  quasi-stationary 
refraction  effects  observed  on  UHF  radar  transmissions  during 
a  three-year  study  at  Millstone  Hill,  Massachusetts. 
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GEOMAGNETIC  LATITUDE 


Figure  33.  The  Latitude  Dependence  of  Ionospheric  Electron 
Content  is  shown,  for  Data  Recorded  at  Uppsala 
(Lietinger  et  al,  Ref.  51). 
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differential-Doppler  and  refraction  error 
perturbations .  The  arrows  mark  the  location 
of  Doppler  excursions  caused  by  gradients 
in  the  high-latitude  ionosphere  (Evans  and 
Wand  (1981),  Ref.  31).  represents  differen¬ 
tial  phase,  and  AEL  represents  elevation  change 
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The  above  discussion  indicates  that  the  trough  is  predomi¬ 
nantly  a  night-time  phenomenon.  However,  Ahmed  and  Sagolyr. 

(Ref.  1)  and  Ahmed  et  al  (Ref.  2)  reported  from  topside 
observations  (1969-1972)  using  ISIS  I  and  INJl^N  V  satellites 
that  in  the  560-3600  km  range,  i.e.  extending  to  altitudes 
well  above  the  peak,  the  trough  had  a  strong  likelihood  of 
being  present  at  any  time:  its  occurrence  frequency  was  95% 
at  night,  60%  near  dawn  and  dusk,  and  48%  near  local  noon. 

They  also  found  negligible  seasonal  variation  in  trough 
location  at  a  given  local  time,  except  near  sunrise. 

(iii)  Equatorward  Displacement  of  Auroral  Phenomena 

It  is  noted  in  Figure  3,  from  Evans  et  al  (Ref.  28) ,  that 
the  Invariant  South  Geomagnetic  Pole  is  displaced  in  the 
direction  of  Australia  by  some  16  deg  of  geographic  latitude 
from  the  South  Geographic  Pole.  As  a  result,  at  the  longi¬ 
tude  of  Adelaide  (138. 6E)  an  invariant  latitude  of  60  S 
corresponds  to  a  geographic  latitude  of  approximately  48  S, 
i.e.  the  mean  trough  position  would  coincide  approximately 
with  the  southern  tip  of  New  Zealand.  Other  invariant 
(southern)  latitudes  (with  geographic  latitudes  indicated 
in  brackets)  are: 

Hobart,  53.6  (42.9);  Macquarie  Island,  64.0  (54.3); 

Mawson,  70.0  (67.4).  (The  co-ordinates  of  many  other  southern 
stations  are  listed  by  Dudeney  and  Piggott  (Ref.  24). 

It  is  also  noted  that  at  a  longitude  of  approximately  142  E, 
the  .South  Geomagnetic  Dip  Pole  is  displaced  23  deg  of  geo¬ 
graphic  latitude  from  the  South  Geographic  Pole  in  the  direc¬ 
tion  of  Australia. 

Coupled  with  the  displacement  of  the  south  magnetic  poles 
towards  Australia  is  the  occurrence  above  the  Southern  Ocean 
of  higher  latitude  ionospheric  phenomena  than  might  otherwise 
be  expected.  In  particular,  the  mean  position  of  the  mid¬ 
latitude  trough  occurs  as  far  north  as  Hobart  in  magnetically- 
disturbed  conditions  (Kp  »  5).  The  trough  commonly  occurs 
at  higher  latitudes  where  there  is,  unfortunately,  a  lack  of 
ionospheric  data. 

6.3  Formation  of  the  Mid-latitude  Trough 

The  mid-latitude  trough  corresponds  to  a  narrow  region  of  very  low 
electron  density.  It  is  situated  just  equatorward  of  the  nocturnal 
auroral  oval  (Figure  32)  within  which  auroral  activity  is  largely 
confined,  and  which  arises  from  diffuse  auroral  particle  precipitation. 

Like  the  auroral  oval,  the  position  of  the  trough  for  a  particular 
degree  of  magnetic  activity  is  "fixed"  in  space  relative  to  the 
:un  while  the  earth  rotates  beneath  it. 

The  trough  is  primarily  a  feature  of  the  night-time  ionosphere: 
it  can  be  very  pronounced  in  winter  and  autumn,  and  less  evident 
in  summer  and  spring. 

The  electron  density  depression  in  the  trough  is  generally  deepest 
near  local  midnight  suggesting  that  the  processes  which  are  acting 
to  destroy  the  ionospheric  plasma  should  be  strongest  at  this  time 
(Schunk  et  al  (Ref.  71)). 
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The  poleward  edge  of  the  trough  is  located  =3  deg  equatorward  of 
the  equatorial  edge  of  the  auroral  oval  (e.g.  Sagolyn  et  al 
(Ref.  66)). 

At  this  geomagnetic  latitude  occurs  the  plasma-pause  which  coincides 
with  the  boundary  between  the  inner  geomagnetic  field  co-rotating 
with  the  eartl  and  the  maqnetospheric  field  resulting  from  the 
geomagnetic  aid  interplanetary  magnetic  fields.  These  are  indicated 
in  Figures  35  and  36,  from  Dudeney  (Ref.  22),  and  Figure  37,  from 
Knecht  (Ref.  46) .  The  envelope  of  the  magnetospheric  field  is 
blown  by  the  solar  wind  to  point  away  from  the  sun. 

The  plasmapause  coincides  with  an  edjrupt  step  in  the  concentration 
of  charged  particles  in  the  atmosphere.  At  lower  latitudes  than 
the  plasmapause  the  atmospheric  plasma  is  relatively  cool  and 
dense;  at  higher  latitudes  the  plasma  is  hot  (i.e.  energetic) 
and  tenuous.  In  the  immediate  vicinity  of  the  plasmapause  is  the 
F  region  mid-latitude  trough,  a  relatively  stable  region  where  the 
concentration  of  ionisation  is  greatly  reduced.  To  date,  no  explana¬ 
tion  proposed  for  the  mid-latitude  trough  is  fully  consistent  with 
the  observational  evidence. 

The  fact  that  the  poleward  wall  of  the  mid-latitude  trough  lies  just 
equatorward  of  the  boundary  of  diffuse  auroral  precipitation  suggests 
that  the  poleward  wall  is  caused  by  ionisation  produced  by  the  soft 
component  of  auroral  particles,  namely  <1  keV  electrons,  precipitated 
from  the  magnetosphere.  The  ionisation  is  produced  at  altitudes 
around  100-150  km,  and  moved  upwards  along  geomagnetic  field  lines 
by  the  effects  of  auroral  heating.  Thus,  the  observed  magnetic 
field-aligned  wall  of  ionisation  is  formed  corresponding  to  a  steep 
horizontal  ionisation  gradient  (as  suggested  in  Figure  36) . 

Statistically,  the  centre  of  the  mid-latitude  trough  lies  near  the 
ionospheric  projection  of  the  geomagnetic  field  line  from  the  equa¬ 
torial  magnetospheric  plasmapause,  sometimes  called  the  magneto¬ 
pause  (Kohnlein  and  Raitt  (Ref.  49)  as  suggested  in  Figure  36. 
However,  at  any  instamt,  the  two  do  not  coincide  i"  a  consistent 
way:  Grebowsky  et  al  (Ref.  36)  explained  this  lack  of  coincidence 
as  follows. 

The  trough  is  coupled  to  the  plasmapause  in  ways  not  fully  under¬ 
stood.  A  very  complex  dynamic  interaction,  in  fact,  occurs.  This 
involves  (i)  the  redistribution  of  plasma  under  the  action  of 
electric  fields  of  magnetospheric  origin  and  thermospheric  winds, 

(ii)  plasma  production  due  to  particle  precipitation  and  solar 
U.V.  radiation,  as  well  as  (iii)  plasma  loss  through  charge  exchange 
and  recombination  processes  (Dudeney  et  al  (Ref.  23)). 

Some  further  detail  on  the  formation  of  the  mid-latitude  trough 
is  given  in  Appendix  3. 
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Figure  35.  Magnetic  Field  Configuration  Near  the  Earth 
(after  Dudeney  (1981),  Ref.  22). 
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Region  of  Reduced  Ionisation  (in  terms  of  elec¬ 
trons/cm*)  emd  its  Location  Relative  to  Halley 
Station  (from  Dudeney  (1981)  ,  Ref.  22)  . 
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7.  DOPPLER  ELECTRON  CONTENT  MEASUREMENTS  -  LOCATING 
THE  TROUGH  IN  REAL  TIME 


The  differential  Doppler  technique  may  be  used  in  measuring  the  spatial 
(and  temporal)  variations  in  total  electron  content  in  the  ionosphere 
between  an  orbiting  satellite  (e.g.  a  NNSS  satellite)  and  a  ground 
station.  The  change  in  phase  difference  between  coherent  150  and  400  MHz 
transmissions  is  measured  at  the  ground  station,  and  is  directly  inter- 
pretcd>le  as  the  change  in  electron  content  along  the  ray  path  from 
satellite  to  ground  station  as  the  satellite  moves  in  its  orbit.  (It 
is  noted  that  phase  changes  are  measured,  and  not  absolute  values  of 
phase  path  length  along  the  ray  path;  i.e.  phase  changes  are  measured 
relative  to  a  minimum  phase,  usually  occurring  at  closest  approach  of 
the  satellite  to  the  receiving  station.  Thus  there  is  always  an  unknown 
constant  =  2Trm  (where  m  is  an  integer)  to  be  added  to  the  measured  phase 
difference  to  obtain  cin  absolute  value  of  phase.) 

Guier  and  Weiffenbach  (Ref.  37)  described  the  Differential  Doppler  tech¬ 
nique  in  proposing  its  use  in  a  satellite  Doppler  navigation  system. 

The  teclmique  is  sumnarised  e.g.  by  Essex  and  Day  (Ref.  27). 

The  differential  Doppler  technique  involves  the  measurement  of  the  differ¬ 
ence  in  phase  path  length  of  radio  signals  received  from  a  satellite  at 
two  harmonically  related  frequencies  f  and  f  * ,  where  f ’  =  mf (m  >  1) . 

Phase  path  expressed  as  a  number  of  wavelengths  may  be  written: 

f 

P  =  -  /  nds 
^  o 

where  n  is  the  refractive  index  of  the  medium  at  a  specified  frequency, 
c  is  the  free  space  speed  of  light,  and 

s  is  the  geometric  path  length  from  satellite  to  receiver. 

Signals  received  at  the  two  frequencies  exhibit  Doppler  phase  path  shifts , 
i.e.  rates  of  change  of  phase  path  length,  of  dP/dt  and  dP'/dt  respec¬ 
tively,  caused  chiefly  by  the  satellite  motion. 


Changes  in  phase  path  length  arising  from  the  line  of  sight  motion  of 
the  satellite  (and  also  phase  delays  introduced  by  the  troposphere  which 
are  frequency  independent)  are  removed  in  calculating  the  differential 
Doppler  frequency  defined  to  be: 


F 


dP 

dt 


dP' 

dt 


(Specifically,  F  is  the  difference  between  the  time  rates  of  change  of 
the  two  phase  paths,  measured  in  a  length  unit  equal  to  the  shorter 
wavelength. ) 

Provided  that  both  frequencies  greatly  exceed  foF2  (the  F2  layer  critical 
frequency) ,  the  refractive  index  of  the  ionosphere  at  a  specified  frequency 
f  has  the  form: 

to  sufficient  accuracy  for  f  >  100  mHz 

(^)'’ 

is  the  ionospheric  plasma  frequency, 

N  being  the  local  electron  concentration,  and 
e  and  m'  the  electron  charge  and  mass  respectively. 
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In  particular,  approximately: 

n  =  1  -  a  jj- 

where  a  is  approximately  4.5  x  10“’  MKS  units. 

It  follows  that  the  differential  Doppler  frequency  (Hz)  is: 

*'=(“-  ^  ^  ^Nsec4>dh 

where  i>  is  the  local  zenith  angle,  and 
hg  the  satellite  altitude. 

Hence : 


/  Fdt  =  /  Nseoljldh  +  c 

It  follows  that  by  integrating  with  respect  to  time  the  differential 
Doppler  frequency  represented  by  P,  it  is  possible  in  principle  to 
determine  the  total  electron  content  (corresponding  to  the  integral  on 
the  right  hand  side  of  the  above  equation)  along  the  ray  path  up  to  the 
satellite  provided  that  the  constant  c  of  integration  can  be  determined. 

Essex  and  Day  (Ref.  27)  refer  to  several  methods  for  determining  c, 
including  the  use  of  local  ionosonde  data  or  combining  differential 
Doppler  data  recorded  at  stations  separated  by  several  hundred  )cm  or  the 
use  of  transmissions  from  a  geostationary  satellite.  However,  in  the 
present  context,  spatial  changes  in  total  electron  content  corresponding 
to  the  presence  of  a  mid-latitude  trough  or  an  atmospheric  disturbance 
can  be  determined  by  differential  Doppler  observations  at  different  times 
(denoted  1  and  2)  at  a  single  station.  In  particular,  the  phase  (number 
of  "cycles")  change  in  the  differential  Doppler  signal  during  a  specified 
time  interval  is  denoted: 


/Padt  -  /Pidt  =  (m  -  /  (Nzsecip^  -  N[Secl^i)dh 

where  the  integral  on  the  right  hand  side  corresponds  to  the  total  electron 
content  change  observed  when  the  satellite  was  in  position  2  compared 
with  position  1. 

Figure  38  shows  exeunples  of  differential  Doppler  frequency  (and  phase) 
versus  latitude  graphs  obtained  by  Essex  and  Day  (Ref.  27) .  Fluctuations 
on  the  frequency  versus  latitude  curves  correspond  to  ionospheric  irregu¬ 
larities. 

Lyon  (Ref.  55)  also  shows  graphs  (Figure  39)  of  electron  content  (along 
a  ray  path)  versus  time,  recorded  over  the  approximately  10  min  intervals 
during  which  NNSS  satellites  are  visible  for  an  earth  observer.  The 
comparatively  smooth  variations  of  electron  content  in  his  graphs  indicate 
the  presence  of  peaics  and  troughs  of  ionisation,  similar  to  those  attribu¬ 
table  to  large  travelling  ionospheric  disturbances  or  the  mid-latitude 
trough.  They  suggest  the  suitability  of  differential  Doppler  observations 
to  record  time  variations  in  electron  content  expected  for  the  mid-latitude 
trough,  as  proposed  in  the  present  project. 

Variations  in  total  ionospheric  electron  content  from  receiver  to  satellite 
are  the  major  cause  of  the  differences  in  phase  variations  observed  at 
two  frequencies,  as  discussed  above.  These  are  the  only  variations  rele¬ 
vant  to  the  present  discussion  which  deals  with  gross  changes  in  electron 
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Figure  38.  Examples  are  shown  of  differential  Doppler 

freauencv  and  phase  versus  latitude  for  throe 
passes  of  the  NNSS  satellites  during  the  period 
29  February  to  1  March,  1976.  Sampling  fret^uency 
IS  2  Hz.  The  times  indicated  are  the  tines  of 
closest  approach  of  the  satellite  (after  Essex 
and  Day  (Ref.  27)). 
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'  e  Tfi-i  fi—1  ft t\ tilde  trouQ^ij  this  effsct  corresponds  to 

density  arising  from  the  Clynch  and  Renfro 

the  second  order  dependence  of  n  on  f,  i-e.  n  i/r  .  y 
(Ref.  19),  however,  mention  other  minor  variations  in  pnase.  namely. 

(a)  a  third  order  term  which  involves  the  geomagnetic  field  ^rough 
the  electron  gyrofreguency  (electron  cyclotron 

this  effect  is  inappreciable  at  frequencies  remote  from  the  gyr 
frequency,  and 

(b)  a  fourth  order  term  which  becomes  appreciable  at  low  elevation  angles. 

CCIR  (Ref.  16)  remarks  that  "maps  of  total  l“ts 

needed  for  predicting  ionospheric  effects  on  the  satellite  Ear^  li^s. 
Measurements  of  total  electron  content  have  been  made  at  many  stations 

S  th.  „«d.y  roftldn  “ 

from  low  orbiting  and  geostationary  satellites  .  Amay  .  - 

»<.  S  pS..„. 

mid- latitude  ionosphere.  Experiments  of  the  ,.ind  describe  P 

re^rt  would  provide  additional  TEC  data  on  ^e  ionosphere  at  high 
invariant  latitudes  (Mendillo  and  Klobuchar  (Ref.  58)). 


Measured  Slant  Content  for  NNSS  Satellite  Pass  at  La  Konge, 
Saskatchewan  on  18  February,  1981,  0341  OT 


32  34  36  36  40 

Tiae.miMitcs 


Measured  slant  Content  and  Content  Fluctuations  for  NNSS  Satellite 
Pass  at  London,  Ontario  on  21  June,  1981,  0234  UT 


Figure  39.  From  Lyon  (Ref.  55). 
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APPENDIX  I 
INVARIANT  LATITUDE 


The  invariant  geomagnetic  latitude  A  is  defined: 

A  =  cos“^(L“'’) 

and  is  based  on  the  pcurameter  L,  defined  by  Mcllwain  (Ref.  56)  as  follows 
(Filz  et  al  (Ref.  33)). 

If  the  geomagnetic  field  were  a  pure  dipole,  L  would  be  equal  to  Rq  which 
is  the  maximum  distance  from  the  centre  of  the  geomagnetic  field  to  the 
farthest  point  on  the  magnetic  field  line  around  which  trapped  charged 
particles  would  be  constrained  to  move. 

It  cem  be  shown  that  in  a  pure  dipole  field,  the  following  relationship 
would  hold: 


where  is  the  magnitude  of  the  magnetic  field  strength  at  the  mirror 
point  of  the  charged  particle  and  is  an  adiabatic  invariant 
of  the  particle  motion, 

M  is  the  magnetic  moment  of  the  dipole, 
f  indicates  a  function,  and 

1  is  the  longitudinal  invariant  which  has  no  clear  physical  or 
geometrical  meaning. 

I  is  defined  as  the  integral  with  respect  to  distance  (from  one  mirror 
point  to  its  conjugate)  of  the  term: 

[‘  -  ©]' 

where  B  is  the  magnitude  of  the  magnetic  field  strength. 

In  the  geomagnetic  field,  the  relationship  (Al.l)  does  not  hold.  Mcllwain 
(1961)  replaced  R^  by  L  which  he  defined  as  follows: 


where  the  units  of  L  are  earth  radii  (Rj.)  ,  and 
M  =  0.311653  (in  units  of  gauss-Rg) . 

Because  the  geomagnetic  field  is  nearly  a  dipole,  L  =  Rq,  as  seen  in 
Figure  40  from  Knecht  (Ref.  46). 

The  locus  of  all  points  with  a  given  L  value  is  called  an  L-shell. 

Invariant  latitude  A  is  a  convenient  geomagnetic  co-ordinate  w)u.ch  gives 
a  good  approximation  to  the  actual  geomagnetic  field,  eund  a  good  indicator 
of  the  behaviour  of  charged  particles  in  the  upper  atmosphere. 

e.g.  For  an  invariant  latitude: 

A  =  cos”*(L  **)  =  60  deg, 
it  is  seen  that: 


i .  e. 


L  =  4 


75 


ERL-0449-SD 


It  follows  that  the  geomagnetic  field  lines  which  emerge  from  the  earth's 
surface  at  an  inveuriant  latitude  of  60  deg,  extend,  very  closely,  to  four 
earth's  radii  from  the  earth's  centre. 

Some  papers  e.g.  Mendillo  and  Chako  (Ref.  57)  employ  corrected  geomagnetic 
co-ordinates.  These  are  defined  (Knecht  (Ref.  46))  by  amalogy  with  geo¬ 
graphic  co-ordinates  emd  represent  the  best  earth-centred  dipole  approxi¬ 
mation  to  the  geomagnetic  field  with  its  axis  tilted  11.5  deg  to  the 
rotation  axis  so  as  to  intersect  the  (currently  approximate)  position 
of  the  geomagnetic  north  pole  at  geographic  co-ordinates  78.5  N,  291  E. 


M. 
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APPENDIX  2 

THE  DOPPLER  MEASUREMENT  OF  RANGE 


The  following  description  of  the  Doppler  technique  in  measuring  the 
receiver-to-satellite  range,  and  hence  in  determining  the  receiver 
position,  is  similar  to  that  given  by  Stansall  (Ref.  76) . 

By  receiving  the  navigation  message,  the  Transit  system  user  learns  the 
position  of  the  satellite  as  a  function  of  time.  Thus,  to  obtain  a 
position  fix  he  must  relate  his  position  to  the  known  satellite  orbit. 
This  relationship  is  established  by  measuring  the  Doppler  shift,  which 
is  a  unique  function  of  the  observer's  position  and  motion  relative  to 
the  known  satellite  orbit. 

Figure  41  illustrates  the  Doppler  measurement  technique.  The  frequency 
fjj  being  received  from  the  satellite  consists  of  the  frequency  being 
transmitted  fq.  plus  a  Doppler  frequency  shift  of  up  to  ±8  kHz  due  to 
relative  motion  between  the  satellite  and  the  receiver.  Note  that  the 
transmitted  frequency  is  offset  low  by  about  80  parts  per  million 
(32  kHz  at  400  MHz)  to  prevent  f^  from  crossing  400  MHz. 

The  navigation  receiver  is  equipped  with  a  stable  reference  oscillator 
from  which  a  400  MHz  ground  reference  frequency  fg  is  derived.  Oscil¬ 
lator  stability  must  be  adequate  to  assume  a  constant  frequency  through¬ 
out  the  satellite  pass.  As  shown  by  Figure  41,  the  navigation  receiver 
forms  the  difference  frequency  fg  -  fj^,  and  each  Doppler  measurement 
is  a  count  of  the  number  of  difference  frequency  cycles  occurring 
between  time  marks  received  from  the  satellite.  Because  every  message 
bit  effectively  represents  another  time  mark,  the  Doppler  counting 
intervals  are  formed  with  respect  to  the  particular  message  format. 

For  example,  each  line  of  the  message  lasts  about  4.6  seconds,  and  the 
conrvonly  used  Doppler  count  interval  of  23  seconds  is  formed  by  starting 
a  new  count  at  the  end  of  every  fifth  line. 

Each  Doppler  count  is  composed  of  two  parts:  the  count  of  a  constant 
difference  frequency  fg  -  frp,  minus  the  count  of  the  number  of  Doppler 
cycles  received  during  that  time  interval.  It  is  the  Doppler  cycle 
count  which  is  physically  meaningful.  The  count  of  the  difference 
frequency  is  an  additive  constant  which  is  eliminated  by  the  position 
fix  calculation. 

In  Figure  41  it  is  emphasised  that  the  distance  between  the  satellite 
and  the  observer  changes  throughout  the  satellite  pass.  It  is  tliis 
change,  in  fact,  which  causes  the  Doppler  frequency  shift.  As  the 
satellite  moves  closer,  more  cycles  per  second  must  be  received  than 
were  transmitted  to  account  for  the  shrinking  number  of  wavelengths 
along  the  propagation  path.  For  each  wavelength  the  satellite  moves 
closer,  one  additional  cycle  must  be  received.  Therefore,  the  Doppler 
frequency  count  is  a  direct  measure  of  the  change  in  distance  between 
the  receiver  and  the  satellite  over  the  Doppler  count  interval.  In 
other  words,  the  Doppler  count  is  a  geometric  measure  of  the  range 
difference  between  the  observer  and  the  satellite  at  two  points  in 
space,  accurately  defined  by  the  navigation  message. 

This  is  a  very  sensitive  measure  because  each  count  represents  one  wave¬ 
length,  which  at  400  MHz  is  only  0.75  metre. 
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The  equation  defining  the  Doppler  count  of  fg  -  fjj  is  the  integral  of 
this  difference  frequency  over  the  time  interval  between  receipt  of 
time  marks  from  the  satellite.  For  example: 

Ni  -  J  (fg  -  f{,)at 

ti  +  Ri/C 

Note  that  ti  +  Ri/C  is  the  time  of  receipt  of  the  satellite  time  mark 
that  was  transmitted  at  time  ti .  The  signal  is  received  after  propa¬ 
gating  over  distance  Rj  at  the  velocity  of  light  C. 

The  cibove  equation  represents  the  actual  measurement  made  by  the  satel¬ 
lite  receiver,  but  it  is  helpful  to  expand  this  equation  into  two  parts; 

N,  =  fgdt  -  f f^dt 

ti  +  Ri/C  ti  +  Ri/C 

Because  the  first  integral  in  this  equation  is  of  a  constant  frequency 
fQ,  it  is  easy  to  integrate,  but  the  second  integral  is  of  the  changing 
frequency  fj^.  However,  the  second  integral  represents  the  number  of 
cycles  received  between  the  times  of  receipt  of  two  timing  marks. 

By  a  "conservation  of  cycles"  argument,  this  quantity  must  equal  identi¬ 
cally  the  number  of  cycles  transmitted  during  the  time  interval  between 
transmission  of  these  time  marks.  Using  this  identity,  this  equation 
can  be  written: 

N  =  f^,dt  -  f^dt 

ti  +  Ri/C  ti 

Because  the  frequencies  f^  and  fy  are  assumed  constant  during  a  satel¬ 
lite  pass,  the  integrals  in  this  equation  become  trivial,  resulting  in: 

Ni  =  fQ  (t2  -  ti)  +  i(R2  -  Ri)j  -  fT(t2  -  ti) 

Rearranging  the  terms  in  this  equation  gives; 

tf 

Ni  =  (fQ  -  f^)  (t2  -  ti)  +  .^(R2  -  Ri) 

This  final  equation  clearly  shows  the  two  parts  of  the  Doppler  count. 

First  is  the  constant  difference  frequency  multiplied  by  a  time  interval 
defined  by  the  satellite  clock.  Second  is  the  direct  measure  of  slant 
range  change  measured  in  wavelengths  of  the  ground  reference  wavelength 
C/fg.  It  happens  that  the  wavelength  of  Eq  is  the  proper  scale  factor 
because  received  time  marks  are  used  to  start  and  stop  the  Doppler  counts. 

If  a  ground  clock  is  used  to  control  the  count  intervals ,  the  wavelength 
of  frj.  would  become  the  appropriate  scale  factor. 

A  usable  satellite  pass  will  be  above  the  horizon  between  10  and  18  minutes, 
which  determines  the  number  of  Doppler  counts  acquired.  Typically  20  to 
40  counts  will  be  collected  by  modern  equipment.  The  Doppler  counts  and 
the  satellite  navigation  message  are  passed  to  a  small  digita’  computer 
for  processing.  For  simplicity,  a  stationary  receiver  is  assumed  as 
shown  in  Figure  41  in  order  to  establish  the  basic  position  fix  concept. 
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APPENDIX  3 

MID-LATITUDE  TROUGH  FORMATION  -  DETAILED  CONSIDERATIONS 


This  appendix  is  intended  to  discuss  in  some  detail  the  formation  of  the 
mid-latitude  trough  as  an  extension  of  its  description  given  in  Section 

4.4. 


Ctoservations  of  the  trough  have  mainly  been  carried  out  by  means  of 
ionospheric  topside  investigations.  During  August  1972  ,  Grebowsky 
et  al  (Ref.  36)  employed  Explorer  45  and  ISIS  2  topside  data,  obtained 
with  a  magnetic  ion  mass  spectrometer  in  studying  the  movement  of  ions, 
particularly  H'*’,  along  geomagnetic  field  lines  associated  with  the  mid¬ 
latitude  trough.  They  proposed  an  explanation  for  the  short-term  differ¬ 
ences  between  the  positions  of  the  trough  and  plasmapause  which  involved 
the  movement  of  ions  at  supersonic  speeds  for  part  of  their  paths  along 
geomagnetic  field  lines.  In  this  case,  ion  densities  would  be  enhanced 
near  the  equatorial  plane  by  means  of  the  (observed)  upward  transport 
of  light  ions  from  the  ionosphere;  nevertheless,  to  agree  with  observa¬ 
tion  the  ionospheric  ion  densities  on  the  same  field  lines  must  remain 
undisturbed.  In  this  way  the  plasmapause  boundary  can  be  formed  at  L 
co-ordinates  (see  Appendix  1)  which,  in  the  equatorial  plane,  are  higher 
than  the  ionosphere.  In  supersonic  ion  flows,  ion  acoustic  waves  cannot 
propagate  "upstream".  Information  about  ion  density  increases  occurring 
at  altitudes  above  where  the  ion  upward  flow  becomes  supersonic  cannot 
therefore  be  propagated  down  to  ionospheric  heights.  In  order  to  account 
at  all  local  times  for  the  separation  between  trough  and  plasmapause, 
Grebowsky  et  al  (Ref.  36)  concluded  that  supersonic  flows  must  persist 
on  the  outer  field  lines  of  the  plasmasphere  for  at  least  one  day  following 
a  storm  depletion  of  ionisation  in  the  trough.  They  concluded  that  upward 
accelerated  H'*’  flows,  possibly  to  supersonic  speeds  during  the  refilling 
of  the  magnetic  flux  tubes  in  the  outer  plasmasphere,  could  produce  an 
equatorial  plasmapause  with  magnetic  field  lines  which  mapped  into  the 
ionosphere  at  latitudes  poleward  of  the  H'*'  density  increase. 

Schunk,  Banks  and  Raitt  (Ref.  70)  concluded  that  there  is  no  single 
cause  for  the  observed  troughs,  but  that  at  various  times,  different 
mechanisms  act  together  to  create  electron  density  depressions  of  substan¬ 
tial  magnitude.  Some  of  the  proposed  mechanisms  are  considered  briefly 
below. 

One  propxjsal  (Nishida  (Ref.  61))  was  that  the  F  region  depletion  could 
be  a  consequence  of  the  escape  of  ionisation  along  open  geomagnetic 
field  lines. 

Knudsen  (Ref.  47)  proposed  that  a  factor  in  the  formation  of  the  trough 
was  a  slow  plasma  convection  across  the  pxjle.  This  enhanced  the  rate 
of  the  ion-atom  interchange  reaction  of  0'*’  with  N^*  (discussed  below) 
resulting  in  a  depletion  of  ionisation  (Schunk  et  al  (Ref.  70) ) . 

The  plasma  convection  pattern  is  illustrated  in  Figure  42  from  Lockwood 
(Ref.  53).  Thermospheric  neutral  winds,  transporting  vibrationally 
excited  N2  molecules  (denoted  N2*)  equatorwards  from  the  auroral  zone, 
would  also  enhance  the  conversion  of  O'*'  to  NO'*'  (Schunk  and  Banks  (Ref. 

69)).  These  factors,  coupled  with  the  normal  lack  of  night-time  main¬ 
tenance  of  the  F  region,  through  recombination,  could  contribute 
appreciably  to  the  observed  reduction  by  an  order  of  magnitude  in  the 
peak  electron  densities  observed  in  the  trough  (Banks  et  al  (Ref.  8); 
Schunk,  Raitt  and  Banks  (Ref.  71)). 
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Figure  42.  The  diagram  is  derived  from  Lockwood  (Ref.  53)  and 

shows  invariant  latitude:  local  time  plots  of  the  model 
plasma  convection  pattern  above  200  km  (from  Knudsen 
(Ref.  48)  )  and  the  mean  positions  of  the  walls  of 
the  trough  in  the  F2-peak  density  at  low  Kp  (between 
the  solid  lines)  obtained  from  regression  equations 
by  Halcrow  and  Nisbet  (Ref.  39) . 
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Through  various  processes,  discussed  by  Schunk  and  Banks  (Ref.  69) ,  N2* 
is  produced  in  the  nocturnal  auroral  oval.  (These  processes  include 
the  excitation  of  N2  by  precipitating  electrons,  thermal  electrons 
and  excited  atomic  oxygen  (O) ,  and  the  reaction  of  N  with  NO.)  The 
presence  of  N2  substantially  increases  the  rate  at  which  O’*"  is  lost 
via  the  reaction: 

0+  +  N2*  NO'*'  +  N  .  (A3.1) 

Owing  to  the  action  of  the  thermospheric  winds  (having  equatorwards 
components)  and  diffusion,  Nj*  does  not  remain  at  its  point  of  origin, 
but  is  transported  equatorwards  to  regions  outside  the  auroral  oval 
where  it  can  act  to  reduce  substantially  the  F-region  electron  density, 
following  the  above  reaction  and  the  one  below  (Spiro  et  al  (Ref.  75)): 

NO'*'  +  e  -*•  N  +  O 

As  Na*  leaves  the  precipitation  region,  there  occurs  a  gradual  loss  of 
vibrational  quanta  and  a  consequent  gradual  decrease  in  the  enhanced  N2 
vibrational  temperature  (T^) .  This  decrease  is  primarily  due  to  a 
"quenching"  reaction  of  excited  O  with  N2,  producing  N2*.  In  the  vici¬ 
nity  of  the  F2  region  peak  electron  density  around  200  km  altitude,  where 
the  quenching  and  diffusion  time  constants  are  equal,  the  characteristic 
time  (T^)  for  the  decay  of  T^  is  probably  determined  by  the  N2*  -  O 
quenching  rate.  (The  loss  of  N2*  by  the  reaction  (A3.1)  above  is  negligible 
compared  with  atomic  oxygen  quenching.)  It  follows  that  is  ~  5000  s 
at  altitudes  above  200  km,  and  if  the  equatorward  component  of  the  neutral 
wind  is  -  100  m/s,  the  vibrating  N2  molecules  travel  approximately  4  deg 
of  latitude  before  quenching  is  likely.  As  the  N2*  ions  drift  over  this 
distance,  which  is  about  the  average  width  of  a  trough,  T^  (previously 
enhanced  in  the  auroral  region)  will  cause  an  enhancement  ir  the  O'*'  loss 
rate  in  reaction  (A3.1)  above.  This  can  lead  to  a  substantial  reduction  in 
the  F  region  ion  density. 

It  is  noted  that  the  cibove  process  involving  N2*  undoubtedly  acts  jointly 
with  other  mechanisms.  Thus  it  can  be  expected  that  within  the  trough, 
transport  processes  will  involve  convection  produced  by  electric  fields 
(discussed  below)  and  thermospheric  winds  with  equatorward  components. 

In  addition,  within  the  trough  the  O'*"  loss  rate  will  be  increased  owing 
to  changes  in  the  reaction  rate  of  equation  (A3.1)  associated  with  electric 
fields  and  N2  vibrational  excitation,  as  well  as  increased  N2*  densities 
arising  from  enhanced  auroral  bombardment  in  the  auroral  oval  (Schunk 
and  Banks  (Ref.  69)). 

The  mid-latitude  trough  can  arise  through  a  depletion  of  ionisation 
involving  both  ionic  reactions,  such  as  those  mentioned  above,  and  plasma 
transport. 

Knudsen  (Ref.  47)  emphasised  a  different  aspect  of  the  nature  of  high- 
latitude  plasma  convection.  In  his  view,  the  low  ion  concentrations 
observed  within  the  trough  resulted  not  from  the  enhanced  recombination 
in  this  region,  but  from  the  longer  time  available  for  recombination  to 
oedur.  , 

Spiro  et  al  (Ref.  75)  also  proposed  that  the  low  ion  concentrations 
observed  within  the  trough  resulted  from  the  long  time  spent  by  plasma 
in  passing  through  this  region,  rather  than  enhanced  recombination  in 
the  region.  Plasma  flow  in  the  high  latitude  ionosphere  is  characterised 
by  a  two-cell  convection  pattern  (seen  in  Figures  42  and  43)  that  results 
from  a  viscous-like  interaction  of  the  wind  with  the  earth's  magneto- 
spheric  boundary  (Axford  and  Hines  (Ref.  7)).  In  pjurticular,  this 
interaction  establishes  a  magnetospheric  polarisation  electric  field. 
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(Another  electric  field  can  possibly  arise  in  the  magnetosphere  from 
the  bulk  motion  with  velocity  V  of  the  solar  wind  past  the  earth  of 
magnetic  induction  g,  creating  a  "V  x  g”  electric  field.  A  further 
electric  field  can  arise  from  a  merging  of  magnetic  field  lines  in 
the  magnetospheric  tail,  creating  an  electric  field  through  the  motion 
of  the  magnetic  field  lines  (Dungey  (Ref.  25)).  By  considering  the 
motion  of  plasma  in  the  electric  fields  and  the  geomagnetic  field,  the 
"tear-drop"  shape  of  the  equatorial  cross-section  of  the  plasmapause 
(Figure  44)  has  been  explained  (e.g.  Kavanagh  et  al  (Ref.  43)). 

In  summary,  Knudsen  (Ref.  47)  appears  to  have  been  the  first  to  offer 
an  explcuiation  for  the  formation  of  an  ionospheric  trough  on  the  night- 
side  of  the  earth.  According  to  this  model,  magnetic  flux  tubes  are 
convected  from  the  dayside  to  the  nightside  across  the  polar  cap  as 
a  result  of  the  in^iressed  magnetospheric  electric  field.  These  tubes 
return  to  the  dayside  by  subsequently  drifting  westwards  or  eastwards. 
Tubes  drifting  westwards  must  travel  faster  than  earth  rotation  if  they 
are  to  re-emerge  in  sunlit  regions  and  some  fail  to  do  so,  allowing 
the  ionisation  in  them  to  fall  to  very  low  levels  as  it  recombines. 

In  tliis  model,  the  trough  then  is  the  locus  of  flux  tubes  which  spend 
a  long  period  in  darkness  owing  to  the  competing  effects  of  the  convec¬ 
tion  electric  field  motion  and  earth  rotation.  This  model  has  been 
refined  further  through  numerical  studies  (Knudsen  et  al  (Ref.  48)) 
and  a  more  realistic  electric  field  convection  pattern  has  been  proposed 
by  Spiro  et  al  (Ref.  75)  based  upon  Atmospheric  Explorer  C  data. 

According  to  the  latter,  no  special  convection  electric  field  pattern 
is  required,  only  that  the  electric  fields  fall  reasonably  smoothly 
as  geomagnetic  latitude  decreases.  There  then  is  always  a  region  on 
the  nightside  where  convection  and  corotation  exactly  balance,  creating 
a  stagnation  line.  Incoherent  scatter  studies  of  the  convection  electric 
field  have  been  made  from  Millstone  Hill  (Evans  et  al  (Ref.  30);  Evans 
et  al  (Ref.  29))  which  appear  to  support  the  explanation  offered  by 
Spiro  et  al  (Ref.  75)  for  the  night-time  trough  as  it  is  a  region  in 
which  the  drift  velocities  appear  to  be  extremely  low. 

Equatorweird  of  approximately  55  deg  invariant  latitude,  the  ionsophere 
is  effectively  shielded  from  the  influence  of  this  magnetospheric 
electric  field.  At  these  lower  latitudes,  a  corotation  electric  field 
dominates.  This  field  arises  through  the  rotation  of  the  dipole  geomag¬ 
netic  field  which  produces  a  radial  electric  field  falling  off  as  1/R^ , 

R  being  the  earth's  radius  (Hones  and  Bergeson  (Ref.  42)).  Plasma  flux 
tubes  within  this  region  corotate  eastward  with  the  earth  under  the 
control  of  the  corotation  electric  field. 

Spiro  et  al's  (Ref,  75)  observations  of  total  ion  concentration  involved 
the  use  of  a  retarding  potential  analyser  on  the  Atmospheric  Explorer  C 
satellite.  Their  data  showed  that  in  the  premidnight  local  time  sector, 
the  plasma  in  the  poleward  portion  of  the  trough  drifts  generally  west¬ 
ward  whie  the  low-latitude  portion  rotates  eastward.  They  explained 
th'-se  observations  in  terms  of  a  plasma  drift  velocity  expressed  as; 

E  X  B 

» 

where  B  represents  the  earth's  magnetic  field  strength,  and 
E  the  local  electric  field. 

The  reversal  of  plasma  drift  velocity  would  arise  through  the  corotating 
electric  field  (effective  at  lower  latitudes)  having  a  direction  opposite 
to  the  magnetospheric  electric  field.  The  long  time  spent  by  plasma  in 
the  trough  region  arises  as  follows.  In  the  premidnight  sector,  plasma 
convection  in  the  poleward  portion  of  the  trough  is  mainly  westward 
along  the  length  of  the  trough  with  only  a  small  equatorward  component. 

In  contrast  in  the  equatorward  portion,  convection  is  mainly  eastward 
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Figure  44. 


AVERACC  CHANGE  M  POSITION  OF  MAIN  TPOUGH 
AT  HALLEY  BAY  (L»t  21 
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The  diagram  is  from  Rodger  and  Pinnock  (Ref.  63) . 
Bars  show  the  mean  velocity  of  main  trough  in  the 
local  hour  immediately  preceding  the  time  when  the 
poleward  edge  of  the  main  trough  is  established 
over  Halley  Bay  (L  =  4.2),  the  locus  of  wtiich  is 
shown  by  the  broken  line.  The  radial  axis  indicates 
the  L  shell  positions.  The  locus  of  the  'tear-drop' 
model  of  the  plasmapause  for  quiet  magnetic  activity 
is  indicated  by  the  continuous  line. 
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along  the  trough  length.  In  summary,  Spiro  et  al  (Ref.  75)  indicated 
that  low  ion  densities  within  the  trough  arise  where  corotation  and 
convective  electric  fields  oppose  each  other,  giving  long  periods  of 
time  in  which  the  plasma  decays. 

Particle  precipitation  and  the  subsequent  redistribution  of  the  resulting 
plasma  by  the  action  of  the  convection  electric  field  were  shown  to  be 
of  primary  infxjrtance  in  the  formation  of  the  poleward  edge  of  the  mid¬ 
latitude  trough  (Pike  et  al  (Ref.  62);  Knudsen  (Ref.  47);  Spiro  et  al 
(Ref.  75);  Soika  et  al  (Ref.  74)).  The  location  of  the  plasmapause  is 
mainly  controlled  by  the  magnitude  of  the  convection  electric  field 
(Lemaire  (Ref.  50))  but  often  there  is  also  significant  particle  preci¬ 
pitation  on  or  near  the  geomagnetic  field  lines  which  thread  the  plasma- 
pause  where  favourable  conditions  exist  for  gyro-resonance  (Gendrin 
(Ref.  34)).  Rodger  and  Pinnock  (Ref.  63)  indicate  that,  for  these 
reasons  a  close  relationship  might  be  expected  between  the  positions  of 
the  trough  edge  and  the  plasmapause.  They  observed  similarities  in  their 
movements  and  positions  on  a  time-scale  of  one  hour  or  more  during  periods 
of  low  geomagnetic  activity.  However,  there  were  dissimilarities  during 
periods  of  high  geomagnetic  activity,  as  was  also  observed  by  Grebowsky  et  al 
(Ref.  36) .  The  similarities  probably  reflect  the  importance  of  the  con¬ 
vection  electric  field  and  particle  precipitation  on  the  formation  of 
the  plasmapause  and  the  poleward  edge  of  the  trough. 

Schunk  et  al  (Ref.  72)  have  combined  a  simple  plasma  convection  model 
with  an  ionospheric-atmospheric  composition  model  incorporating  ideas 
of  the  kind  outlined  above  in  this  section,  in  order  to  simulate  high- 
latitude  ionospheric  behaviour  which  is  found  to  be  in  reasonable  agree¬ 
ment  with  observation. 
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APPENDIX  4 

COMPUTATION  OF  THE  LOCATION,  HEIGHT  AND  lATITUDIRAL  EXTENT  OF  THE 
MID-LATITUDE  TROUGH 


This  appendix  describes  the  results  of  preliminary  coipputations  of  the 
location,  height  and  latitudinal  extent  of  the  mid-latitude  trough,  and 
briefly  indicates  the  reasoning  {following  that  in  Section  4.5)  which 
led  to  the  computations  which  were  performed  by  Mr.  G.D.  Slade. 

A  modified  version  of  the  DRCS  RADFIN  program  was  employed.  The  lati¬ 
tudinal  extent  of  the  trough  was  determined  using  a  subroutine  developed 
by  B.W.  Halcrow  and  J.S.  Nisbet  based  on  their  paper  (Ref.  39). 

The  height  of  the  F2  layer  in  the  region  where  the  trough  exists  was 
determined  by  modifying  the  CCIR  Worldwide  Numerical  Map  Data  (Ref.  15). 
From  the  discussion  in  Section  4.4,  the  height  hj^  of  the  contour  of  maximum 
F2  region  electron  density  within  the  trough  was  assumed  to  be  200  km 
above  the  height  h^,  of  an  equal  electron  density  contour  in  the  CCIR 
(i.e.  "normal")  data; 

i.e.  hj,  =  hg  +  200  km 

As  indicated  in  Section  4.5,  the  ratios  of  the  heights  of  other  electron 
density  contours  within  the  trough  were  taken  to  be  the  same  as  the  ratios 
for  the  corresponding  electron  density  contours  derived  from  Ref.  15. 

(In  this  context  the  "height"  was  taken  to  be  the  altitude  of  a  contour 
within  the  trough  minus  the  altitude  of  the  contour  remote  from  the  trough.) 

An  example  of  preliminary  computations  made  on  this  basis  is  shown  in 
Figure  45. 
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Figure  45.  A  computed  raid- latitude  trough  is  shown  superimposed  on  contours  of  constant 
electron  density  (having  plasma  frequencies  (MHz)  as  indicated)  derived  from 
data  in  Kef.  15  for  the  longitude  of  Adelaide  at  11.00  UT  in  July  and  zero 


ERL-0449-SD 


88  - 


REFERENCES 


No.  Author  Title 

1  Ahmed,  M.  and  "Topside  Ionospheric  Trough  Morphology  at 

Sagolyn,  R.C.  Mid-  and  High-Latitudes", 

Compilation  of  papers  presented  by  the 
Space  Physics  Division  at  the  Ionospheric 
Effects  Symposium,  Air  Force  Geophysics  Lab., 
Hanscon  AFB,  Massachusetts, 

Report  AFGL-TR-78-0080,  April  1978. 

2  Ahmed,  M. ,  "Topside  Ionospheric  Trough  Morphology: 

Sagolyn,  R.C.,  Occurrence  Frequency  and  Diurnal,  Seasonal 

Wildmcin,  P.J.L.  and  and  Altitude  Variations", 

Burke,  W.J.  J.Geophys.Res. ,  489,  1979. 

3  Alfven,  H,  "Some  Properties  of  Magnetospheric  Neutral 

Surfaces" , 

J.Geophys.Res.,  73,  13,  1968. 

4  Alouette  I  "Alouette  1  Ionospheric  Data  -  Interpolated 

N(h)", 

1-6,  Department  of  Communications,  Defense 
Research  Telecommunications,  Ottawa,  Canada, 
1968. 

5  Aluoette  II  "Alouette  II  Ionospheric  Data  -  Interpolated, 

N(h)", 

1-2,  Department  of  Conanunications ,  Defense 
Research  Telecommunications,  Ottawa,  Canada, 
1970. 

6  Amayenc,  P. ,  Bertin,  F,  "Sur  I'Evolution  Latitudinal  de  Contenu 
and  Papet-Lepine,  J.  Electronique  de  1 ’ Ionosphere" , 

Ann. da  Geophys.,  27_,  345,  1971. 

7  Axford,  W.I.  and  "A  Unifying  Theory  of  High-Latitude 

Hines,  C.O.  Phenomena  and  Geomagnetic  Storms", 

Can.J.Phys.,  39,  1433,  1961. 

8  Banks,  P.M.,  "NO'*'  and  0+  in  the  High-Latitude  F-region" , 

Schunk,  R.W.  ana  Geopiiy s. Rev.  Lett. ,  239,  1974. 

Raitt,  W.J.  . 

9  Belrose,  J.S.  "The  Lower  Ionospheric  Regions", 

Chapter  3  of  "Physics  of  the  Earth's  Upper 
Atmosphere”,  ed.  Hines,  C.O.,  Paghis,  J. 
Kartz,  T.R.  and  Fejer,  J.A.,  1965 

10  Bowman,  G.G.  "Directional  Characteristics  of  lonosonde 

Interference  Patterns  from  the  Filchner  Ice 
Shelf", 

J. Atmos. Ter r .Phys. ,  30 ,  1115,  1968. 

11  Bowman,  G.G.  "Ionization  Troughs  Below  the  F2-layer 

M.  imum" , 

Planet. Space  Sci.,  jj,  777,  1969. 


89 


ERL-0449-SD 


Author 

Buchau,  J.,  Pike,  C.P. 
and  Wong ,  M . 


Burge,  J.D., 

King,  J,W.  and 
Slater,  A.J. 

C.C.I.R, (International 
Radio  Consultative 
Committee) 

C.C.I.R. (International 
Radio  Consultative 
Committee) 


C.C.I.R.  (International 
Radio  Consultative 
Committee) 

Chako,  C.C.  and 
Mendillo,  M. 

Chan,  K.L.  and 
Colin,  L. 


Clynch,  J.R-  and 
Renfro,  B.A. 


COMSAT  General 
Corporation 

Conkright,  R.O.  and 
Brophy,  H,J. 


Dudeney,  J.R. 


Dudeney,  J.R., 
Jarvis,  M. J. , 
Kressman,  R.I., 
Pinnock,  M. , 
Rodger,  A.S»  and 
Wright,  K.H. 


Title 

"Detailed  Specification  of  the  Arctic 
Ionosphere  and  an  Application  of  Three- 
Dimensional  Raytracing" , 

Air  Force  Cambridge  Research  Laboratories, 
Report  No.  AFCRL-TR-7 3-0726,  November  1973. 

"Mapping  of  ^F2  by  Means  of  Topside 
Sounder  Satellites", 

TeleconBiunication  Journal,  40,  356,  1973. 

Reconmendations  and  Reports  of  the  C.C.I.R., 
"Propagation  in  Ionized  Media" , 
vol.  VI,  XIVth  Plenary  Assembly  (Kyoto) , 
Geneva,  1978. 

Report  340  concerning  the  C.C.I.R.  Atlas 
of  Ionospheric  Characteristics  (Oslo  1966) , 
International  Telecommunication  Union, 
Geneva,  1967. 

15th  Plenary  Assembly  (Document  6/1004-e) , 

Geneva,  1982. 


Report  No.  AFGL-TR-7B-0092(I) ,  Air  Force 
Geophysics  Labs.,  Hanscomb,  .Massachusetts,  1967. 

"Global  Electron  Density  Distributions  from 
Topside  Soundings" , 

Proc.  IEEE,  990,  1969. 

"Evaluation  of  Ionospheric  Range  Error  Model" , 
Proc.  Third  International  Geodetic  Symposium 
on  Satellite  Doppler  positioning,  Vol.  1; 
Sponsor;  Dense  Mapping  Agency  and  National 
Ocean  Survey?  Host:  Physical  Science  Lab.  of 
New  Mexico  State  Uni.,  Las  Cruces,  497, 

February,  1982. 

"MARISAT  System  Description", 
prepared  by  C.G.C.  as  Manager  of  the  MARISAT 
Joint  Venture,  March  1977. 

Catalog  of  Ionosphere  Vertical  Soundings 
Data,  Report  UAG-85,  world  Data  Center  A 
for  Solar-Terrestrial  physics,  U.S.  Dept, 
of  Commerce  ,  National  Oceanic  and  Atmos¬ 
pheric  Administration  Environmental  Data 
and  Information  Service,  Boulder,  Colorado, 

U.S. A.,  80303,  July  1982. 

"The  Ionosphere  -  a  View  from  the  Pole", 

New  Scientist,  September  17,  1981. 

"Ionospheric  Troughs  in  Antarctica", 

Nature,  295,  307,  1982. 


ERL-0449-SD 


90 


Title 


No .  Author 


24 


25 


26 


27 


28 


29 


30 


31 


32 


Dudeney,  J.R.  and  “Antarctic  Ionospheric  Research", 

Piggott,  W.R.  Antarctic  Research  Series  (American 

Geophysical  Union) ,  29,  200,  1978. 


Dungey,  J.W.  "Interplanetary  Magnetic  Field  and  the 

Auroral  Zones” , 

Phys. Rev, Letters,  47,  1961. 


ERNO  Raumfahrttechnik  "SARSAT,  Ein  Rettungssystem  fur  Schiffe 
GMBH  und  Flugzeuge", 

Report  80-054,  Bremen,  West  Germany, 
October  1980.  (English  translation: 
NASA  Technical  Memo.,  NASA  TM-76391, 
"SARSAT:  A  Rescue  System  for  Ships  and 
Airplanes" . 


Essex,  E.A.  and 
Day,  T.W. 


Evans,  J.E. 
Newkirk,  L.L.  and 
KcCormac,  B.w. 


"The  Differential  Doppler  Technique  and  its 
Use  in  Measuring  the  Spatial  and  Temporal 
Variation  of  Irregularities  in  the 
Ionosphere" , 

Monitor  -  Proc.  of  the  IREE  Aust. ,  104 ,  1978. 

"North  Polar,  South  Polar,  World  Maps  and 
Tables  of  Invariemt  Magnetic  Co-ordinates 
for  Six  Altitudes:  0,  100,  300,  600,  1000 
and  3000  km" , 

Final  Report  DASA  2347  (N-9A-69-1)  to  the 
Defense  Atomic  Support  Agency,  Washington  D.C. 
from  Lockheed  Palo  Alto  Research  Laboratory', 
October,  1969. 


Evans,  J.V.,  Holt,  J.M.,  "Millstone  Hill  Incoherent  Scatter  Observa- 
Oliver,  W.L.  and  tions  of  Auroral  Convection  over 

Ward,  R.H.  60  deg  Invariant  Latitude  75  deg, 

2.  Initial  Results", 

J.Geophys.Res. ,  85,  41,  1980. 

Evans,  J,V,,  Holt,  J.M.,  "Millstone  Hill  Incoherent  Scatter  Observa- 
and  Wand,  R.H.  tions  of  Auroral  Convection  over 

60  deg  Invariant  Latitude  ^75  deg, 

I,  Observing  and  Data  Reduction  Procedures", 

J. Geophys.Res.,  8^,  7059,  1979. 


Evans,  J.V.  and  "Anomalous  Ionospheric  Refraction  Associated 

Wand,  R.H.  with  the  Auroral  Zone", 

article  in  "Effect  of  the  Ionosphere  on 
Radiowave  System" (editor-in-chief  J.M.  Goodman) 
based  on  Ionospheric  Effects  Symposium 
at  Old  Town,  Alexandria,  Vancouver,  sponsored 
by  Naval  Res.  Lab.  Office  of  Naval  Res.  and 
Air  Force  Geophys.  Lab.,  April  1981. 


Feinblum  D.A.  and  "HILION  -  A  Model  of  the  High-Latitude 

Horan,  R.J.  Ionospheric  f2  Layer  and  Statistics  of 

Regular  Ionospheric  Effects  at  Fort  Churchill, 
1968" , 

Joint  Radar  Propagation  Study  of  Bell 
Laboratories,  Western  Electric  and  Lincoln 
Laboratory,  Massachusetts  Institute  of  Tech¬ 
nology,  for  U.S.  Army  Safeguard  System 
Command  under  Contract  DAHC60-71-C-00G5  ur  !  for 
Advanced  Ballistic  Missile  Defence  Agency 
under  Contract  F19628-7 3-C-OOC2 ,  1973. 


91 


ERL-0A49-SD 


No .  Author  Title 

33  Filz,  R.C.,  Katz,  L. ,  "Corpuscular  Radiation", 

Kuck,  G.A.,  Shea,  M.A.  A  revision  of  Chapter  17,  "Handbook  of 
and  Smart,  D.F.  Geophysics  and  Space  Environments", 

Air  Force  Cambridge  Research  Labs., 
AFCRL-68-0666,  Air  Force  Surveys  in 
Geophysics,  No.  208,  December  1968. 

34  Gendrin,  R.  '  Is  the  Plasmapause  a  Preferential 

Region  for  Proton  Precipitation?", 
Annls.Geophys. ,  127,  1975. 

35  George,  P.L.  “HF  Skywave  Aspects  of  Trunk-Area 

Communications  ESM", 

Department  of  Defence,  Defence  Science 
and  Technology  Organisation,  Electronics 
Resecirch  Laboratory,  Defence  Research 
Centre,  Salisbury,  Tech.  Memo, 
ERL-0234-TM  (unclassified),  April  1982. 


36 

Grebowsky,  J.M., 

"Ionospheric  and  Magnetospheric  Plasma- 

Hoffman,  J.H.  and 

pauses"  , 

Maynard,  N.C. 

Planet. Space  Sci.,  26,  651,  1978. 

37 

Guier,  W.H.  and 

"A  Satellite  Doppler  Navigation  System" 

Weiffenbach,  G.C. 

Proc.  I.R.E.,  507,  1960. 

38 

Halcrow,  B.W. 

"F2  Peak  Electron  Densities  in  the  Main 

Trough  Region  of  the  Ionosphere" , 

Report  No.  PSU-IRL-IR-55 ,  Contract: 

NASA  NGL  39-009-003,  Ionosphere  Research 
Lab.,  Pennsylvania  State  Uni.,  May  1976. 

39  Halcrow,  B.W.  and  "A  Model  of  F2  Peak  Electron  Densities  in 

Nisbet,  J.S.  the  Main  Trough  Region  of  the  Ionosphere", 

Radio  Sci.,  1^,  815,  1977. 

40  -  "Hcmdbook  for  Radiotelephone  Ship  Station 

Operators" , 

ISBN  0-642-01431-0  (Australian  Government 
Printing  Service,  Canberra),  1981. 

41  Heaps,  M.G.  "Accounting  for  Ionospheric  Variability 

and  Irregularity  in  High  Frequency  Direction 
Finding" , 

article  in  "Effect  of  the  Ionosphere  on 
Radiowave  Systems"  (editor-in-chief 
J.M.  Goodman)  based  on  Ionospheric  Effects 
Symposium  at  Old  Town,  Alexandria, 

Vcmcouver,  sponsored  by  Naval  Res.  Lab. 
Office  of  Naval  Res.  and  Air  Force 
Geophys.  Lab.,  April  1981. 

42  Hones,  E.W.  cind  "Electric  Field  Generated  by  a  Rotating 

Bergeson,  J.E.  Magnetized  Sphere", 

J. Geophys. Res. ,  70,  4951,  1965. 

43  Kavanagh,  L.D.  "Plasma  Flow  in  the  Magnetosphere", 

Freeman,  S.W.  and  J. Geophys. Re s. ,  73 ,  5511,  1968. 

Chen,  A.J. 


Author 


92 


Title 


ERL-0449-SP 

tio. 


44 


45 


46 


47 


48 


49 


50 


51 


52 


53 


54 


55 

56 


King,  J.W.  "The  Determination  of  FqF2  and  hn,F2  from 

Satellite-borne  Probe  Data" , 
Telecommunication  Journal,  40,  364,  1973. 

Klobuchar,  J.A.  “Total  Electron  Content  Studies  of  the 

Ionosphere" , 

Report  TR-73-0098,  Air  Force  Cambridge 
Resecirch  Labs.,  Bedford,  Massachusetts, 
U.S.A.,  1973. 


Knecht,  D.J.  "The  Geomagnetic  Field  (A  Revision  of 

Chapter  11,  Hamdbook  of  Geophysics  and 
Space  Environments", 

Air  Force  Surveys  in  Geophysics  No.  246, 
AFCRL-72-0570,  September  1972. 


Knudsen,  W.C.  “Magnetospheric  Convection  and  the  High- 

Latitude  F2  Ionosphere" , 

J.Geophys.Res. ,  79,  1046,  1974. 


Knudsen,  W.C.,  "Numerical  Model  of  the  Convecting  f2 

Bcinks,  P.M.  Ionosphere  at  High  Latitudes", 

Winningham,  J.D.  and  J.Geophys.Res,,  82 ,  4784,  1977. 
Klumpeir,  D.M. 


Kohnlein,  W.  and 
Raitt,  W.J. 


Lemaire,  J. 


Lietinger,  R., 
Hedberg,  A.  and 
Tanskanen ,  P . 


"Position  of  the  Mid-Latitude  Trough  on  the 
Topside  Ionosphere  as  Deduced  from  Esro  4 
Observations" , 

Planet. Space  Sci.,  25,  600,  1977. 

"The  Mechanisms  of  Formation  of  the  Plasma- 
pause"  , 

Annls.Geophys. ,  31,  175,  1975. 

"The  Electron  Content  of  the  Ionosphere  and 
the  Southern  Boundary  of  Diffuse  Aurora", 
J. Atmos. Terr .Phys. ,  44,  369,  1982. 


Liszka,  L. 


Lockwood,  M. 


Lockwood,  M. 


"The  Htgh-Latitude  Trough  in  Ionospheric 
Electron  Content", 

J. Atmos. Terr. Phys. ,  29^,  1243,  1967. 

"The  Bottoraside  Mid-Latitude  Ionospheric 
Trough" , 

J. Atmos. Terr .Phys. ,  42 ,  605,  1980. 

"A  Simple  Model  of  the  Effects  of  the  Mid- 
Latitude  Total  Ion  Trough  in  the  Bottomside 
F  layer  on  HF  Radiowave  Propagation" , 

Radio  Science,  16,  3,  385,  1981. 


Lyon,  G.F. 


"Ionospheric  Effects  of  the  SARSAT  System", 
Canad, Aeronaut. Space  Journal,  27 ,  327,  1981. 


Mcllwain,  C.E. 


Mendillo,  M.  and 
Chako,  C.C. 


"Co-ordinate  for  Mapping  the  Distribution  of 
Magnetically  Trapped  Particles" , 
J.Geophys.Res.,  66,  11,  3681,  1961. 

"The  Baselevel  Ionospheric  Trough" , 
J.Geophys.Res.,  82^,  5129,  1977. 


57 


No. 


Author 


93 


Title 


ERL-0449-SD 


58 


59 


60 


61 


62 


63 


64 


65 


66 


67 


68 


69 


70 


Mendillo,  M.  and 
Klobuchcu:,  J.A. 


Huldrew,  D.B. 


Mundo,  C.,  Tami,  L. 
and  Larsen,  G. 


Nishida,  A. 


Pike,  C.P., 
Whalen,  J.A.  eind 
Buchau,  J. 

Rodger ,  A . S .  and 
P innock,  M. 


Rodger,  A.S.  and 
Pinnock,  M. 


Rush,  C.M. , 
Rosich,  R.K. , 
Brooks,  C.B.  eind 
Leise,  D.L. 


Rycroft,  M.J.  and 
Burnell,  S.J. 


Rycroft,  H.J.  eind 
Thomas,  J.O. 


Sagolyn,  R. 
Bredesen,  S.  and 
Wlldman,  P. 


Schunk,  R.W.  cuid 
Banks,  P.M. 


Schunk,  R.W., 
Banks,  P.M.  ^md 
Raitt,  W.J. 


"Investigations  of  the  Ionospheric  F-region 
Using  Multi-station  Total  Electron  Content 
Observations" , 

J.Geophys.Res.  ,  80^,  643,  1975. 

"F-layer  Ionization  Troughs  Deduced  from 
Alouette  Data" , 

J.Geophys.Res.,  70 ,  IS,  3635,  1965. 

"Final  Report  Program  Plan  for  Search  and 
Rescue  Electronics  Alerting  and  Locating 
System", 

DOT-TSC-OST-73-42,  February  1974. 

"Average  Structure  of  the  Storm-time  Change 
in  the  Polar  Topside  Ionosphere  at  Sunspot 
Minimum" , 

J.Geophys.Res.,  7^,  6051,  1967. 

"A  12-hour  Case  Study  of  Auroral  Phenomena 
in  the  Midnight  Sector" , 

J.Geophys.Res.,  82^,  3547,  1977. 

"The  Variability  and  Predictability  of  the 
Main  Ionospheric  Trough”, 
from  C.S.  Deehr  and  J.A.  Holder  (eds.), 
"Exploration  of  the  Polar  Upper  Atmosphere" 

(D.  Reidel  Publ.  Co.),  463,  1980. 

"Movements  of  the  Mid-Latitude  Ionospheric 
Trough" , 

J. Atmos. Terr. Phy s.  ,  985,  1982. 

"A  Simplified  Model  of  the  High  Latitude 
Ionosphere  for  Telecommunications 
Applications" , 

NTIA  Report  82-94,  U.S.  Dept,  of  Commerce, 
National  Telecommunications  emd  Information 
Service,  Boulder,  Colorado,  Janueury  1982. 

"Statistical  Analysis  and  Movements  of  the 
Ionospheric  Trough  and  the  Plasmapause" , 
J.Geophys.Res.,  75 ,  5600,  1970. 

"The  Magnetospheric  Plasmapause  and  the 
Electron  Density  Trough  at  the  Alouette  I 
Orbit” , 

Planet. Space  Sci.,  65,  1970. 

Article  in  "An  Empirical  Model  of  the  Polar 
Ionosphere",  ed.  T.J.  Elkins, 

Air  Force  Surveys  in  Geophysics  No.  267, 
AFCRI/-TR-73-0331,  .lay  1973. 

"Auroral  Nj  Vibrational  Excitement  of  the 
Electron  Density  Trough" , 

Geophys. Res. Letters,  6,  239,  1975. 

"Effects  of  Electric  Fields  and  Other  Processes 
Upon  the  Night-Time  High-Latitude  F-layer", 
J.Geophys.Res.,  81,  3271,  1976. 


Author 


94 


Title 


ERL-044g-SD 

No. 


71  Schunk,  R.W., 
Raitt,  W.j.  and 
Banks,  P.M. 

72  Schunk,  R.W., 
Raitt,  W.J.  and 
Sojka 


73  Sharp,  G.W. 


74  Soika,  J.J., 

Raitt,  W.J.  and 
Schunk ,  R . W . 


75  Spiro,  R.W. , 
Heelis,  R.A.  and 
Hanson,  W.B. 

76  Stansell,  T.A. 


77  Thomas,  J.O.  and 
Dufour,  S.W. 

78  Tulunay,  V.K. 


79  Tulunay,  Y.K.  and 
Grebowsky,  J.M. 


80  Tulunay,  Y.K.  and 
Sayers,  J. 


81  Wrenn,  G.I,.  and 
Raitt,  W.J. 


Zawalick,  E.J.  and 
Cage,  A.L. 


"Effect  of  Electric  Fields  on  the  Day-time 
High-Latitude  E  and  F  regions", 

J.Geophys.Res. ,  80,  3121,  1975. 

"High-Latitude  Ionospheric  Model:  First 
Step  Towards  a  Predictive  Capability" , 
article  in  "Effect  of  the  Ionosphere  on 
Radiowave  Systems"  (editor-in-chief 
J.M.  Goodman)  based  on  Ionospheric  Effects 
Symposium  at  old  Town,  Alexandria, 

Vancouver,  sponsored  by  Naval  p-^c  T,ab. 

Office  of  Naval  Res.  and  Air  Force 
Geophys.  Lab.,  April  1981. 

"Mid-Latitude  Trough  in  the  Night  Ionosphere" , 
J.Geophys.Res.,  'H,  1345,  1966. 

"Theoretical  Predictions  for  Ion  Composition 
in  the  High-Latitude  Winter  F-region  for 
Solar  Minimum  and  Low  Magnetic  Activity", 
J.Geophys.Res.,  86,  2206,  1981. 

"Ion  Convection  and  the  Formation  of  the 
Mid-Latitude  F  Region  Ionization  Trough", 

83,  A9,  4255,  1978. 

"The  Transit  Navigation  Satellite  System", 
Magnavox  Government  and  Industrial 
Electronics  Co.,  R-5933,  October  1978. 

"Electron  Density  in  the  Whistler  Medium" , 
Nature,  London,  206,  567,  1965. 

"Magnetically  Symmetric  Detection  of  the  Mid- 
Latitude  Electron  Density  Trough  by  Ariel  III 
Satellite", 

J. Atmos. Terr .Phys. ,  34,  1547,  1972. 

"The  Noon  and  Midnight  Mid-Latitude  Trough 
as  seen  by  Ariel  IV" , 

J. Atmos. Terr .Phys. ,  40,  845,  1978. 

"Characteristics  of  the  Mid-Latitude  Trough 
as  Determined  by  the  Electron  Density 
Experiment  of  Ariel  III", 

J. Atmos. Terr .Phys. ,  33,  1737,  1971. 

"In  Situ  Observations  of  Mid-Latitude 
Ionospheric  Phenomena  Associated  with 
the  Plasmapause" , 

Ann.(3eophys. ,  31,  17,  1975. 

"Frequency  Tables  of  the  Geomagnetic  Index 
Kp,  1932-1970", 

J.Geophys.Res.,  76,  7009,  1971. 


82 


ERL-0449-SD 


DISTRIBUTION 


EXTERNAL 

In  Unlt«d  Kingdom 

Counsellori  Defence  Science,  London 

British  Library,  Document  Supply  Centre  (UK) 

Institution  of  Electrical  Engineers 

British  Antarctic  Survey,  Natural  Environ  Res  Council 
(Attention:  Dr  J.R.  Dudeney) 

(Attention:  Dr  W.R.  Plggott) 

University  of  Leicester,  Ionospheric  Physics  Group 
(Attention:  T.B.  Jones) 

Rutherford-Appleton  Laboratory 
(Attention:  Dr  J.W.  Ring) 

University  College  London,  Mullard  Space  3cl  Lab 
(Attention:  Dr  W,J.  Raltt) 

University  of  Southampton,  Dept  of  Physics 
(Attention:  Dr  M,J.  Rycroft) 

University  College  of  Vales,  Physics  Dept 
(Attention:  Prof  L.  Thomas) 

In  United  States  of  America 

Counsellor,  Defence  Science,  Washington 

NASA  Scientific  and  Technical  Information  Office 

Regis  College 

(Attention:  Dr  M.  Ahmed) 

Pennsylvania  State  University,  Ionospheric  Res  Lab 
(Attention:  Dr  R.L.  Dlvany) 

SAMSO 

(Attention:  Dr  B.W.  Halcrow) 

Lockheed  Palo  Alto  Res  Lab 

(Attention:  Dr  W.C.  Rnudsen) 

Boston  University,  Dept  of  Astronomy 
(Attention;  Dr  M.  Mendlllo) 

Pennsylvania  State  University,  Ionospheric  Res  Lab 
(Attention:  Dr  J.S.  Nlsbet) 

Air  Force  Cambridge  Res  Lab,  Ionospheric  Physics  Lab 
(Attention:  Dr  R.C.  Sagolyn) 


Copy  No. 

Cnt  Sht  Only 
1 
2 

3 

4 

5 

6 

7 

8 

9 

Cnt  Sht  Only 
10 

U 

12 

13 

14 

15 

16 


17 


ERL-0449-SD 


University  of  California  at  San  Diego,  Dept  Applied  Phys 

(Attention:  Dr  R.W.  Sc hunk)  18 

Magnavox  Gov  and  Indust  Electronics 

(Attention:  Mr  T.A.  Stansell)  19 

In  Canada 

University  of  Western  Ontario,  Physics  Dept 

(Attention:  Prof  G.F,  Lyon)  20 

In  Getiaany 

University  of  Bonn,  Astronomical  Institute 

(Attention:  Dr  W.  Rohnlein)  21 

In  Turkey 

Technological  University 

(Attention:  Dr  Y.K.  Tulunay)  22 

In  New  Zealand 

University  of  Auckland,  Radio  Res  Centre 

(Attention:  Dr  M.  Lockwood)  23 

In  Australia 

Chief  Defence  Scientist 

Assistant  Chief  Defence  Scientist  (Policy)  24 

Assistant  Chief  Defence  Scientist  (Operations) 

Joint  Intelligence  Organisation  (DSTI) 

(Attention;  HST-AF)  25 

Defence  Signals  Directorate  26 

Director  General,  Joint  Communications  and  Electronics 

(Attention:  DJCE(PP))  27 

Director,  Naval  Intelligence  and  Security 

(Attention:  DDNI)  28 

Director,  Electronic  Warfare  -  Navy  29 

Director,  Military  Intelligence 

(Attention:  SOl(EW))  30 

Director  General,  Operational  Requirements  -  Air  Force 

(Attention:  OREW)  31 

Director,  Air  Force  Intelligence  and  Security 

(Attention;  ISINT)  32 

Director,  Communications  and  Electronics  -  Air  Force 

(Attention:  CERCOM2)  33 


ERL-0449-SD 


Document  Exchange  Centre, 

Defence  Information  Services  Branch  for: 

Microfiche  copvlng  34 

United  Kingdom,  Defence  Research  Information  Centre  35  -  36 

United  States,  Defense  Technical  Information  Center  37  -  48 

Canada,  Director,  Scientific  Information  Services  49 

New  Zealand,  Ministry  of  Defence  50 

National  Library  of  Australia  51 


Director  General,  Army  Development  (NSO) ,  Russell  Offices 


for  ABCA  Standardisation  Officers 

UK  ABCA  representative,  Canberra  52 

US  ABCA  representative,  Canberra  53 

Canada  ABCA  representative,  Canberra  54 

NZ  ABCA  representative,  Canberra  55 

Librarian,  Technical  Reports  Centre,  Defence  Central 
Library,  Campbell  Park  56 

Library,  DSD,  Melbourne  57 

Library,  Materials  Research  Laboratories  58 

Library,  RAN  Research  Laboratory  59 

Library,  Aeronautical  Research  Laboratories  60 

Director,  Industry  Development,  Adelaide  Cnt  Sht  Only 

South  Aust  Inst  of  Tech,  Physics  Dept 

(Attention:  Prof  G.L.  Goodwin)  61  -  62 

University  of  Queensland,  Physics  Dept 

(Attention:  Dr  G.G.  Bowman)  63 

University  of  Adelaide,  Physics  Dept 

(Attention:  Dr  B.H.  Briggs)  64 

(Attention:  Dr  R.A.  Vincent)  65 

LaTrobe  University,  Dlv  of  Theoretical  and  Space  Physics 

(Attention:  Prof  K.D.  Cole)  66 

INTERNAL 

Director,  Electronics  Research  Laboratory  67 

Superintendent,  Electronic  Warfare  Division  68 

Superintendent,  Communications  Division  69 

Principal  Officer,  Radio  Wave  Propagation  Group 


70 


ERL-0449-SD 


Dr  K.J.W.  Lynn,  Radio  Wave  Propagation  Croup 

Library,  Defence  Science  and  Technology  Organisation 
Salisbury 

CD  Documentation  Section 
Spares 


71 

72  -  73 
74  -  75 
76  -  78 


DOCUMENT  CONTROL  DATA  SHEET 


UNCLASSIFIED 


2  ;  SECURITY  CLASSIFICATION 
;  a.  Complete 

;  Document:  Unclassified 


b.  Title  in 

Isolation:  Unclassified 


i  I  c.  Summary  in 

I  j  Isolation:  Unclassified 


3  '  DOWNGRADING /DELIMITING  INSTRUCTIONS 

Limitation  to  be  reviewed  in 
March  1991 


4  I  TITLE  I 

LOCATING  THE  SOURCE  OF  LONG-DISTANCE  RADIO  DISTRESS  SIGNALS  FROM  THE 
SOUTHERN  OCEAN 


5  !  PERSONAL  AUTHOR  (S) 


G.L.  Goodwin 


8  i  8.  1  CORPORATE  AUTHOR  (S) 


Electronics  Research  Laboratory 


8.2  DOCUMENT  SERIES 
and  NUMBER 
Special  Document 
0449 


11  I  IMPRINT  (Publishing  organisation) _ 

Defence  Science  and  Technology 
Organisation  Salisbury 


13  I  RELEASE  LIMITATIONS  (of  the  document) 


Approved  for  Public  Release. 


Security  classification  of  this  page :  | _  UNCLASSIFIED 


6  DOCUMENT  DATE 
March  1988 


TJ  7. 1  TOTAL  NUMBER 

OF  PAGES  94 


7.  2  NUMBER  OF 

REFERENCES  8? 


9  REFERENCE  NUMBERS 


a.  Task ; 


b.  Sponsoring  Agency ; 


10  '  COST  CODE 


720350 


12  '  COMPUTER  PROGRAM  (S) 
(Title  (s)  and  language  (s)) 


Security  classification  of  this  page : 


1  !  DOCUMENT  NUMBERS 
j  AR 

]  Number:  AR-005-365 


Series 

Number:  ERL-0449-SD 


Other 
Numbers : 


Security  classification  of  this  page : 


UNCLASSIFIED 


14  ANNOUNCEMENT  LIMITATIONS  (of  the  inlofinalion  on  thasa  pages) 


No  limitation 


15  DESCRIPTORS  16  COSATI  CODES 

a.  EJC  Thesaurus  Position  (location) 

Terms  Radio  signals  0045C 

Search  and  rescue 
Ionospheric  propagation 


b.  Non  -  Thesaurus 

Terms  NAVSAT 


17  SUMMARY  OR  ABSTRACT 

(if  this  Is  security  classifisd.  the  announcement  of  this  report  enll  be  similarly  classified) 


This  report  summarises  a  feasibility  study  and  strategy  appropriate  to 
search  and  rescue  operations  for  ships  or  aircraft  in  distress  in  the 
Southern  Ocean.  The  use  of  high  frequency  ionospherical ly-propagated 
radio  transmissions  is  justified.  Measurement  of  the  direction  of 
arrival  of  a  radio  signal  facilitates  ray-tracing  through  the  ionosphere 
Characteristics  of  the  F2  region  mid-latitude  ionisation  trough  are 
considered  in  detail.  Because  the  trough  is  displaced  towards  Australia, 
its  position  must  be  known  as  near  as  practicable  in  real  time  for 
sufficiently  accurate  ray-tracing  to  be  performed  quickly  enough  after 
receipt  of  a  distress  signal.  The  differential  Doppler  technique  using 
e.g.  150  and  400  MHz  tranmissions  tram  NAVSAT  (NNSS)  satellites  is 
recommended  in  locating  the  trough  position. 


Security  classification  of  this  pag*:  |  UNCLASSIFIED 


I 


